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Abstract 

Compelling evidences in favor of neutrino masses and mixing obtained 
in the last years in Super-Kamiokande, SNO, KamLAND and other 
neutrino experiments made the physics of massive and mixed neutrinos 
a frontier field of research in particle physics and astrophysics. There 
are many open problems in this new field. In this review we consider 
the problem of the absolute values of neutrino masses, which appar¬ 
ently is the most difficult one from the experimental point of view. We 
discuss the present limits and the future prospects of /3-decay neutrino 
mass measurements and neutrinoless double-/? decay. We consider the 
important problem of the calculation of nuclear matrix elements of neu¬ 
trinoless double-/? decay and discuss the possibility to check the results 
of different model calculations of the nuclear matrix elements through 
their comparison with the experimental data. We discuss the upper 
bound of the total mass of neutrinos that was obtained recently from 
the data of the 2dF Galaxy Redshift Survey and other cosmological 
data and we discuss future prospects of the cosmological measurements 
of the total mass of neutrinos. We discuss also the possibility to obtain 
information on neutrino masses from the observation of the ultra high- 
energy cosmic rays (beyond the GZK cutoff). Finally, we review the 
main aspects of the physics of core-collapse supernovae, the limits on 
the absolute values of neutrino masses from the observation of SN1987A 
neutrinos and the future prospects of supernova neutrino detection. 
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1 Introduction 


Compelling evidences in favor of neutrino oscillations, driven by small neutrino masses 
and neutrino mixing, were obtained in the Super-Kamiokande [1 r2l 1 3l!3lEJ , SNO 001 , 
KamLAND 0 and other atmospheric mum, solar iiaimaim and long-baseline HH! 
neutrino experiments. These findings have brought the physics of massive and mixed 
neutrinos in the front line of the research in particle physics and astrophysics 1 . 

From all the existing terrestrial and astrophysical data it follows that neutrino masses 
are smaller than the masses of other fundamental fermions (lepton and quarks) by many 
orders of magnitude. There is a general consensus that the smallness of neutrino masses 
is due to New Physics beyond the Standard Model. In the most attractive see-saw 
mechanism of neutrino mass generation PI. the smallness of neutrino masses is due 
to the violation of the total lepton number on a scale which is much larger than the 
clectroweak scale. 

There are many open problems in the physics of massive and mixed neutrinos: 

• How many light neutrinos with definite mass exist in nature? 

The minimal number of massive neutrinos is equal to the number of active (flavor) 
neutrinos (three). If, however, sterile neutrinos exist, the number of massive neu¬ 
trinos is larger than three (see Refs. prpnp. The data of all the existing neutrino 
oscillation experiments (solar [T^ll31ll4llT5|IH10 0l£] . atmospheric and 

LSND [21]) require the existence of (at least) four massive neutrinos. LSND is the 
single accelerator experiment in which the transition > P e has been observed. 
The check of the LSND claim is an urgent problem. This will be done by the 
MiniBooNE experiment Pi. which started recently. 

• What is the nature of the neutrinos with definite mass: are they purely neutral 
Majorana particles, or Dirac particles, possessing a conserved total lepton number? 

The answer to this fundamental question can be obtained through the investigation 
of processes in which the total lepton number is not conserved. The most promising 
process is neutrinoless double /3-decay of some even-even nuclei. There are many 
new experiments on the search for neutrinoless double /3-decay now in preparation 
(see Ref. j23|h which will push the experimental sensitivity at a level that is about 
two orders of magnitude better than today’s sensitivity. We will discuss neutrinoless 
double /3-decay in this review. 

• What are the absolute values of the neutrino masses? 

Neutrino oscillations are due to differences of phases which different massive com¬ 
ponents of the initial flavor neutrino states pick up during their evolution. As a 
result, neutrino oscillation experiments allow to obtain information only on neu¬ 
trino mass-squared differences (see Refs. ( 2111231111120 ] )• If is very important that 
neutrino oscillation experiments are sensitive to tiny neutrino mass-squared differ¬ 
ences, because of the possibility to explore very large distances and small energies. 
The measurement of the absolute values of neutrino masses at a level of a few eV 

1 See Ref. TJ, for an extensive bibliography on neutrino physics and astrophysics. 
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is a challenging problem. This review is dedicated to the discussion of this problem 
(see also Ref. ). 

Let us mention also the very important problems of the precise determination of the 
values of the neutrino oscillation parameters and the search for CP violation in neutrino 
oscillations. These problems will be investigated in experiments at the future Super-Beam 
facilities and Neutrino factories (see Refs. We will not discuss them here. 

We will start in Section |21 with a short review of the present status of neutrino os¬ 
cillations. We will consider neutrino oscillations in solar, atmospheric and long-baseline 
neutrino experiments in the framework of mixing of three neutrinos. The importance for 
neutrino mixing of the results of the long-baseline reactor experiments CHOOZ and Palo 
Verde, in which no indication in favor of neutrino oscillations was found, will be stressed. 

In Section 01 we will consider the Mainz mmm and Troitsk pT^lfTTTTj experiments 
on the measurement of the neutrino mass through the detailed investigation of the end¬ 
point part of the /9-spectrum of tritium. We will discuss also the future KATRIN tritium 
experiment [551. 

In Section 01 we briefly review the most recent results of the experiments on the 
measurement of the effect of neutrino masses in pion and tau decays. 

Section 01 is dedicated to neutrinoless double /3-decay. Even though in this review we 
are mainly interested in the possibilities to obtain information about the absolute values 
of the neutrino masses from the investigation of this process, some aspects of the theory 
of the process will be also presented. 

The role of neutrinos in Astrophysics and Cosmology has been under the scrutiny of 
physicists with ever increasing intensity over the last few decades (see Refs. mmm 
mmmmmm)- Actually, the intimate relationship of neutrinos and astrophysics 
goes even further back in time when Bethe and others realized that the inner workings 
of the Sun proceed via the thermonuclear reactions that burn hydrogen into helium and 
release neutrinos. Since then there has been a steady increase in interest and involvement 
in the study of neutrinos in astrophysical environments. Not only the interest in solar 
neutrinos has been particularly intense over the last years, where such epochal events as 
their detection on Earth in a variety of underground experiments have been milestones of 
late 20th century high energy physics, but also the involvement of neutrinos in stellar core 
collapse has been theoretically analyzed and observationally established in the momentous 
detection of neutrinos of SN1987A and, furthermore, the influence and role of neutrinos 
in the cosmic evolution has been a major area of research in contemporary high energy 
physics. 

Seminal work on neutrinos and Cosmology was pursued in the late sixties and early 
seventies when neutrinos appeared as ideal candidates to contribute substantially to the 
matter density of the Universe (see Ref. BSD- In fact, hot dark matter models were 
popular for quite some time as they seemed to render a satisfactory model for structure 
formation. Of course, the interest in neutrinos worked then and still works now both 
ways; from the cosmological arena, neutrinos were welcomed in the new cosmological 
Paradigm but also from the Particle Physics side, Cosmology/Astrophysics was used and 
is used even more so today to constrain and sharpen the still not well known properties 
of neutrinos. 

The main focus of this review is the mass of neutrinos and especially their absolute 
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mass. So we have selected the issues in Cosmology/Astrophysics that have relevant im¬ 
pact on the extraction of information concerning neutrino mass. They are contained in 
Sections El □ and El We start in Subsection EZD by introducing the famous Gerstein- 
Zeldovich upper bound on the total sum of neutrino masses that can contribute to the 
matter density of the Universe mm- Subsection 16.21 is dedicated to an overview of the 
temperature fluctuations of the Cosmic Microwave Background Radiation (CMB) with a 
special attention to the characteristics of the peak structure of the angular power spec¬ 
trum (see Ref. There, the influence of neutrino mass on the anisotropy spectrum 

is explicitly discussed. Although it is shown that this influence is not as significant as 
the role of other cosmological parameters that enter the angular spectrum of temperature 
fluctuations, Subsection Id. 21 is relatively long as compared to the other Subsections in Sec¬ 
tion El because in any analysis of cosmological import the CMB is of pivotal importance. 
Subsection IQ1 is devoted to Galaxy Redshift Surveys. Neutrino mass has a remarkable 
effect on the power spectrum of matter distribution and this effect is observable in the 
large samples of data compiled in present galaxy distribution surveys or to be collected 
in future surveys. The final astrophysical source of information discussed in this review, 
namely Lyman a forests studies, is dealt with in Subsection 16.41 The last Subsection 
in Section El Subsection 16.51 contains the summary of all relevant neutrino mass limits 
obtained in the actual analysis by different groups and by different authors of the astro- 
physical/cosmological sources that have been discussed in the foregoing Subsections, ft 
contains also a brief report on the prospects for neutrino mass in this rapidly changing 
held of Cosmology and Astrophysics. 

Another topic that we cover in our review concerns cosmic rays. A probe of neutrino 
properties could come from the observation of cosmic rays with energies exceeding the 
Greisen-Zapsepin-Kuzmin cutoff mm- A possible explanation could be the so-called 
Z-burst scenario EHE3, where a flux of ultra high energy neutrinos interacts with relic 
cosmological neutrinos, producing cosmic rays through the Z-resonance. The resonance 
condition involves the masses of neutrinos and we review the status of this mechanism in 
Section 0 (see also Ref. [2BJ) • 

In 1987 the observation of neutrinos coming from supernova 1987A in the Large 
Magellanic Cloud marked the beginning of extra solar system neutrino astronomy and 
allowed to get information on the supernova mechanism and neutrino properties (see 
Refs, mm)- particular, the values of the neutrino masses are limited by the lack 
of spread of the observed neutrino signal, which would be caused by energy-dependent 
velocities of sufficiently massive neutrinos [HHlEniEIlEElEHl • In Section El we review the 
classification and rate of supernovae (Section EH}, the current theory of core-collapse su¬ 
pernova dynamics (Section EH}, the observation of SN1987A neutrinos (Section EH}, the 
inferred limits on neutrino masses (Section EH}, and the future prospects for supernova 
neutrino detection (Section EH} ■ 


2 Status of neutrino oscillations 


Strong evidences in favor of neutrino oscillations were obtained recently in Super-Kamiokande 


um\mm 


SNO 


GflHS! 


, KamLAND P and other atmospheric mm, solar mmm 
and long-baseline up neutrino experiments. These findings gave us the first evidence that 
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neutrino masses are different from zero and the fields of neutrinos with definite mass v l 
enter into the standard charged current (CC) and neutral current (NC) 

3a ° = VlLldL 1 ja C = 5^ (2.1) 

l=e,[i,T l=e,n,r 

in the mixed form 

VIL = ^2 U H u iL (l = e, /I, r ), (2.2) 

i 

where U is the unitary mixing matrix. The minimal number of massive neutrinos z/j is 
equal to the number of active (flavor) neutrinos (three). The number of massive neutrinos 
can be larger than three (see Ref. Hi). In this case, in addition to Eq. EH we have 

^sL ^ ^ UsiViL , (2.3) 

i 


where v sR (s = Si, S 2 , • • •) are the fields of sterile neutrinos 2 . 

The most plausible mechanism of neutrino mass generation is the see-saw mechanism 
|IB] • In order to explain this mechanism, let us consider the simplest case of one generation 
and assume that the standard Higgs mechanism with one Higgs doublet generates the 
Dirac mass term 

£ d = -m D v R v L + h.c.. (2.4) 

It is natural to expect that rrin is of the same order of magnitude as the mass of the charged 
lepton or quarks in the same generation. We know, however, from experimental data that 
neutrino masses are much smaller than the masses of charged leptons and quarks. I 11 order 
to “suppress” the neutrino mass let us assume that there is a lepton-number violating 
mechanism beyond the Standard Model which generates the right-handed Majorana mass 
term 3 

£ R 3 =-~ Mr 9 R (v R ) c + h.c. , (2.5) 

with M r ?ud (usually it is assumed that Mr ~ Mgut ~ 10 15 GeV). Here (v R ) c = 
CVr T , where C is the charge conjugation matrix. 

After the diagonalization of the total neutrino mass term, for the light neutrino mass 
we obtain 

1 1 / 

m = -~M r + - (2.6) 

I 11 the case of three generations, the see-saw mechanism leads to a spectrum of masses 
of Majorana particles with three light neutrinos with masses rn t (z = 1, 2, 3) much smaller 
than the quark and charged-lepton masses, and three very heavy masses of the order of 
the scale of violation of the total lepton number (for recent reviews see jflOllfllj i. 

2 The fields v s l do not enter into the standard CC and NC in Eq. O- They could be right-handed 
neutrino fields, SUSY fields, etc.. 

3 Notice that, since v R is a SU(2) singlet and has zero hypercharge, the Majorana mass term £))' is 
allowed by the electroweak gauge symmetries. 


6 





Let us stress that, if the neutrino masses have a standard see-saw origin 4 , neutrinos 
with definite masses are Majorana particles. In some models which implement the see¬ 
saw mechanism (see, for example, Ref. [5H), the neutrino masses naturally satisfy the 
hierarchy 

mi -C m 2 m 3 . (2.7) 

If there is neutrino mixing, the state of a neutrino (active or sterile) with momentum 
p is given by the coherent superposition of the states of neutrinos with definite masses 

, ( 2 - 8 ) 

where \vj) is the state of neutrinos with momentum p, mass ra* and energy 

/ TYl^ 

E i — \jp 2 + mf ~ p+ (p 2 > ml ). (2.9) 

From Eq. (EED it follows that if at time t = 0 a neutrino u a (a = e, //, r) is produced, the 
probability amplitude to find u a i at time t is given by 

A(y a v a ,) = (u a > \e~ lHot \u a ) = U <*'i ^ Ki • (2-10) 

i 

Thus, the probability of the transition v a —> zv in vacuum is given by 

2 

5a'a + UohUY (e~ iAm ^ - l) . (2.11) 

i> 2 



v«) = E K, 

i 


Here Am), = m 2 —m 2 , L ~ t is the distance between the neutrino source and the neutrino 
detector, and E is the neutrino energy. 

In the simplest case of transitions between two types of neutrinos (n /t —> u T or > u e , 
etc.), the index i in Eq. (IXTTD takes only one value i = 2 and for the transition probability 
we obtain the standard expression 

1 — cos Am 2 —;] (a 1 ^ a ), (2-12) 

2 E J 

where Am 2 = m\ — m 2 , |I/ Q 2 | 2 = cos 2 id, \U a i 2 \ 2 = 1 — |I/ a 2 | 2 = sin 2, d (d is the mixing 
angle). For the probability of v a to survive we have 

2 2 1 ) ’ (2 ' 13) 


P(^ Q 


= 1 - P(n Q 


v n ') = 1 


sin 2 2d ( 1 — cos A m 


P(n Q —> v a i) = ^ sin 2 2d 


The expressions (EUH) and (ITT31) describe periodical transitions between two types of 
neutrinos (neutrino oscillations). They are widely used in the analysis of experimental 
data. 

4 In non-standard see-saw models neutrinos could be Dirac particles jS21R)3j. 
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Figure 1: Zenith angle distribution of Super-Kamiokande sub-GeV single ring e-like 
events, //-like events, multi-GeV single ring e-like events, and //-like events + partially 
contained (PC) events. The black histogram shows the Monte Carlo prediction and the 
gray histogram is the best fit for —> u T oscillations with Am 2 = 2.5 x lCC 3 eV 2 and 
sin 2 2i9 = 1.0. Figure taken from Ref. jfiTj . 

2.1 Atmospheric neutrinos 

The first model independent evidence in favor of neutrino oscillations was obtained in the 
atmospheric Super-Kamiokande (S-K) experiment [TpUd] . In this experiment a significant 
zenith angle 9 Z asymmetry of the high energy muon events was observed. At high energies 
the zenith angle 9 Z is determined by the distance L, which neutrinos pass from the 
production region in the atmosphere to the detector. If there are no neutrino oscillations, 
the number of detected multi-GeV (E > 1.3 GeV) electron (muon) events satisfy the 
symmetry relation 

Ni(cos9 z ) = A(-cos6b) (Z = e, /x) (2.14) 

As one can see from Fig. |U the number of multi-GeV electron events observed in the 
S-K experiment is in good agreement with this relation. On the other hand, Fig. |T| shows 
that the multi-GeV muon events observed in the S-K experiment strongly violate the 
relation (ixm Let us define the ratio U/D , where U is the number of up-going muons 
(—1 < cos 6b < —0.2, 500 km < L < 13000 km) and D is the number of down-going 
muons (0.2 < cos 9 Z < 1, 20km < L < 500km). For the multi-GeV events in the S-K 
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Figure 2: Allowed region of the oscillation parameters Am 2 = Am atm and sin 2 26 = 
sin 2 2 , d atm from the analysis of Super-Kamiokande atmospheric neutrino data. Figure 
taken from Ref. [FT] . 


experiment it was found 


(U/D)me as 

(U/D)uc 


0.54 ± 0.04 (stat.) ± 0.01 (syst.), 


(2.15) 


where (U/D) M q is the ratio predicted by Monte Carlo under the assumption that there 
are no neutrino oscillations. If there are no neutrino oscillations, the ratio of ratios 
(U/ D) meas /(U/ D )mc must be equal to one. The S-K value (12.151) differs from one by 
11 a. 

The data of the S-K [HE, (21 and other atmospheric neutrino experiments (SOUDAN 
2 TO] , MACRO jllj ) are well described assuming that two-neutrino oscillations 
v r take place. The allowed region of the neutrino oscillation parameters Am atm and 
sin 2 2 , d atm from the analysis of the S-K data is shown in Fig. [3 At 90% C.L. the oscillation 
parameters are in the ranges 


1.6 x 10~ 3 < Am 2 tm < 3.9 x 10" 3 eV 2 , sin 2 2-d atm > 0.92 . (2.16) 


The best-fit values of the parameters are 


Am 2 tm = 2.5 x 10” 3 eV 2 , sin 2 2tf atm = 1.0 ( X 2 mm = 163.2/170 d.o.f.). (2.17) 


Recently the results of the first long-baseline accelerator experiment K2K have been 
published [ITJ- hi this experiment neutrino oscillations in the atmospheric range of Am 2 
were searched for. Neutrinos mainly from decays of pions, produced by 12 GeV protons 
hitting a beam-clump target at the KEK proton accelerator, were detected by the S-K 
detector at the distance of about 250 km from the source. The average neutrino energy 
is 1.3 GeV. 
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Figure 3: Allowed regions of the oscillation parameters Am 2 = Am^ 2K an d sin 2 26 = 
sin 2 2 '$k 2K obtained from the results of the K2K experiment. The dashed, solid and dot- 
dashed lines are, respectively, 68.4%, 90% and 99% C.L. contours. The best fit point is 
indicated by the star. Figure taken from Ref. lIBI- 


In the K2K experiment there are two near detectors at the distance of about 300 
nr from the beam-dump target: a 1 kt water Cherenkov detector and a fine-grained 
detector. The number and the spectrum of muon neutrinos detected in S-K are compared 
with the expected quantities, calculated on the basis of the results of the near detectors. 
Quasielastic one-ring events + n —> n~ + p are selected for the measurement of the 
energy of the neutrinos. 

The total number of muon events observed in the S-K experiment is 56. The expected 
number of events is 80.l!®' 2 . The observed number of one-ring muon events that was used 
for the calculation of the neutrino spectrum is 29. The expected number of one-ring events 
is 44. 

The regions of the allowed values of the oscillation parameters obtained from a maxi¬ 
mum likelihood two-neutrino oscillation analysis of the K2K data are presented in Fig. 0 
The best-fit values of the parameters are 

sin 2 2'd K2 K = 1, Am| 2K = 2.8 x 1(T 3 eV 2 . (2.18) 

These values are in agreement with the values of the oscillation parameters obtained from 
the analysis of the S-K atmospheric neutrino data (see Eqs. (12.1611 and emu). 

Thus, the K2K experiment confirms the evidence for neutrino oscillations that was 
found in the atmospheric Super-Kamiokande experiment. The K2K data reported in 
Ref. m have been obtained with 4.8 x 10 19 protons on target (POT). The K2K experi¬ 
ment is planned to continue until about 10 2 ° POT will be reached. 
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2.2 Solar neutrinos 


The event rates measured in all solar neutrino experiments are significantly smaller than 
the event rates predicted by the Standard Solar models. The following values were ob¬ 
tained, respectively, for the ratio of the rates observed in the Homestake JT2], GALLEX- 
GNO (Him, SAGE jjBTj and S-K folRTHI experiments and those predicted by the BPOO 0T)J 
Standard Solar Model (SSM): 0.34 ± 0.03, 0.58 ± 0.05, 0.60 ± 0.05, 0.465 ± 0.018. It has 
been known during many years that these data can be explained by transitions of the 
initial solar yds into other neutrinos, which cannot be detected in the radiochemical 
Homestake, SAGE, GALLEX and GNO experiments. The S-K experiment is sensitive 
mainly to yds. 

Recently, strong model independent evidence in favor of the transitions of solar yds 
into y/s and yds has been obtained in the SNO experiment mm- In this experiment 
solar neutrinos are detected via the observation of the following three reactions 5 : 

• The charged-current (CC) reaction 

v e J rd^e~+p + p ) (2.19) 

• The neutral-current (NC) reaction 

v x + d^v x + n + p, (2.20) 

• The elastic-scattering (ES) reaction 

u x T 6 —> y, T e . (2.21) 

The kinetic energy threshold for the detection of electrons in the CC and ES processes 
in the SNO experiment is 5 MeV. The NC process has been detected through the ob¬ 
servation of 7 rays from the capture of neutrons by deuterium. The NC threshold is 2.2 
MeV. Thus, practically only neutrinos from 8 B decay are detected in the experiment 6 * . 

The measurement of the total CC event rate allows to determine the flux of y, on the 

|V o V» 

*7 = - "e))cc < . (2.22) 

where 4> 6 e is the total initial flux of yds and (P(y, —*► y>))cc is the probability of y, 
to survive averaged over the CC cross section and the known initial spectrum of S B 
neutrinos. In the SNO experiment it was found that 

(^ C )sno = (l.76t 6 ; 66 (stat.)t 6 ; 66 (syst.)) x 10 6 cni- 2 S - 1 . (2.23) 

All active neutrinos y>, y, and y- are recorded by the detection of the NC process 
(I TM . Taking into account the universality of neutral currents (see the recent analysis 

5 Here v x stands for any active neutrino. 

6 According to the BPOO SSM gSl, the flux of high energy hep neutrinos is about three orders of 

magnitude smaller than the flux of 8 B neutrinos. Looking at solar neutrino events beyond the 8 B 

spectrum endpoint, the Super-Kamiokande Collaboration found that the flux of hep neutrinos is smaller 
than 7.9 times the BPOO SSM hep flux at 90% C.L. IBBj . 


11 





in Ref. im the total flux of all active neutrinos on the Earth measured in the SNO 
experiment is 


(^ C )sno = Y ( $ ^ C )sno = (5.09^;43(stat.)l^(syst.)) x 10 6 cm 2 s \ (2.24) 

l=e,fi,r 

which is about three times larger than the CC flux in Eq. 

All active neutrinos are detected also via the observation of the ES process (irm 
However, the cross section of the neutral-current z/ /i;T e —> v^ T e scattering is about six 
times smaller than the cross section of the charged-current and neutral-current v e e —> v e e 
scattering. 

The event rate R ES of the ES process (12.2111 can be written as 

R ES = Ke) ■ (2-25) 


Here ( a Vf , e ) is the v e e —> v e e cross section averaged over the initial spectrum of 8 B neu¬ 
trinos and the ES flux <h ES is given by 


$ ES 

V 



(ch/^e) 

W^e) 


$ ES 

i'/i.T ’ 


(2.26) 


where <& Ve and <h ;y r are, respectively, the fluxes of v e and v^ T on the Earth averaged over 
the ES cross section and the initial spectrum of 8 B neutrinos. The ratio of the averaged 
/y /i;T e —> v^ T e and u e e —> v e e cross sections is given by 


7 ’/ 1 &) r\ r' A 

^ ~ 0.154 . 

(^ee) 


In the SNO experiment it was found 

(4 )ES )gNo = (2.39i^23(stat.) ± 0.12(syst.)) x 10 6 cm _2 s _1 . 


(2.27) 


(2.28) 


This value is in good agreement with the value of the ES flux determined in the S-K 
experiment PEHj. In the S-K experiment solar neutrinos are detected via the observation 
of the ES process (irm During 1496 days of running a large number, 22400 ± 800, of 
solar neutrino events with recoil energy above the 5 MeV threshold were recorded (the 
uncertainty is due to the statistical subtraction of background events). From the data of 
the S-K experiment it was obtained 

(4> es )s_k = (2.35 ± 0.02(stat.) ± 0.08(syst.)) x 10 6 cm" 2 s” 1 . (2.29) 


In the S-K experiment also the spectrum of the recoil electrons was measured. No sig¬ 
nificant distortion of the spectrum with respect to the expected one (calculated under 
the assumption that the shape of the spectrum of v e on the Earth is given by the known 
initial 8 B spectrum) was observed. Furthermore, no distortion of the spectrum of the 
electrons produced in the CC process was observed in the SNO experiment. These 

data are compatible with the assumption that the probability of solar neutrinos to survive 
is a constant in the high energy 8 B region. Thus, we have 


djNC _ ^CC _ (f>ES 
v e — ^ v e — ^v e ■ 


(2.30) 
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4> e (10 6 cm" 2 s' 1 ) 

Figure 4: Flux of and v T vs flux of v e in the 8 B energy range deduced from the three 
neutrino reactions in SNO. The diagonal bands show the total 8 B flux as predicted by 
the BPOO SSM [HH! (dashed lines) and that measured with the NC reaction in SNO (solid 
band). The intercepts of these bands with the axes represent the ±lcr errors. The bands 
intersect at the fit values for 0 e = $j, e and T , indicating that the combined flux 

results are consistent with neutrino flavor transformation assuming no distortion in the 
8 B neutrino energy spectrum. Figure taken from Ref. [7j. 


Obviously, the NC flux can be presented in the form 


$ N c = $ 

V V, 


NC 


+ $ 


NC 


(2.31) 


Combining the CC and NC fluxes and using the relation (USD, we can determine now 
the flux In Ref. 0 the ES flux was also taken into account as an additional 

constraint (see Fig. H The resulting flux of and v T on the Earth is 

($v m , t )sno = ( 3 - 41 — 0 I 45 (stat.)^o] 45 (syst.)) x 10 6 cm" 2 s _1 . (2.32) 


Thus, the detection of the solar neutrinos through the simultaneous observation of CC, 
NC and ES processes allowed the SNO collaboration to obtain a direct model independent 
5.3 u evidence of the presence of and u T in the flux of the solar neutrinos on the Earth. 

Before the publication of the first results t6j of the SNO experiment, from the global fit 
of the data of the Homestake [T2] , SAGE |IS], GALLEX |I3j, GNO Tlj and S-K [7T1177] 
experiments several allowed regions in the plane of the two-neutrino oscillation parameters 
Am s 2 ol and tan 2 d so \ had been found: the large mixing angle (LMA), low mass (LOW) 
and small mixing angle (SMA) Mikheev-Smirnov-Wolfenstein (MSW) [73,[74] regions, 
the vacuum oscillation (VAC) region and others (see, for example, Ref. [UJEJt]). The 
situation changed after the publication of the first SNO data [§], which, together with 
the recoil electron spectrum measured in the S-K experiment dunam, disfavored the 
SMA-MSW region (see, for example, Ref. ESI)- The most recent data from the SNO m 
and S-K [UEH] experiments strongly disfavor the SMA-MSW region (see, for example, 
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Figure 5: Left: The ratio of measured to expected v e flux from reactor experiments. 
The shaded region indicates the range of flux predictions corresponding to the 95% C.L. 
LMA region found in a global analysis of the solar neutrino data [ 22 ] • The dotted curve 
corresponds to the best-fit values = 5.5 x 1CT 5 eV 2 and sin 2 2$ so i = 0.83 found in 

Ref. [221 • Right: Upper panel: Expected reactor v e energy spectrum with contributions 
of h geo (antineutrinos emitted by 238 U and 232 Th decays in the Earth) and accidental 
backround. Lower panel: Energy spectrum of the observed prompt events (solid circles 
with error bars), along with the expected no oscillation spectrum (upper histogram, with 
h geo and accidentals shown) and best fit (lower histogram) including neutrino oscillations. 
The shaded band indicates the systematic error in the best-fit spectrum. The vertical 
dashed line corresponds to the analysis threshold at 2.6 MeV. Figures taken from Ref. pj. 

Ref. [77.). All global analyses of the present solar neutrino data favor the LMA-MSW 
region IHUZHHZZUZmimEIIlEIHHlllHn] • The best-fit values of the oscillation parameters in 
the LMA-MSW region found in Ref. J8j are 

Am 2 ol = 5 x 1(T 5 eV 2 , tan 2 d sol = 3.4 x KT 1 , (y% n = 57.0 , 72 d.o.f.) (2.33) 

In the next Subsection we will discuss the recent results of the long-baseline reactor 
experiment KamLAND j2j. The data of this experiment allow to exclude the SMA, 
LOW and VAC regions, leaving the LMA region as the only viable solution of the solar 
neutrino problem. 

2.3 The first results of the KamLAND experiment 

Recently the first results of the KamLAND experiment, started in January 2002, have 
been published [2]. In this experiment electron antineutrinos from many reactors in Japan 
and Korea are detected via the observation of the process 

v e + p —■> e + + n . (2.34) 

The threshold energy of this process is E * h ~ m e + m n — m p = 1.8 MeV. About 80% of 
the v e flux is expected from 26 reactors with distances in the range 138-214 km. The 1 kt 
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Figure 6: Left: KamLAND excluded regions of neutrino oscillation parameters Am 2 = 
^ m KamLAND an d sin 2 20 = sin 2 2d Ka mLAND for the rate analysis and allowed regions for 
the combined rate and energy spectrum analysis at 95% C.L. At the top are the 95% 
C.L. excluded region from CHOOZ jHlj and Palo Verde 1551 experiments, respectively. 
The dark area is the 95% C.L. LMA allowed region obtained in Ref. IS2). The thick dot 
indicates the best fit of the KamLAND data in Eq. (12.361) . Figure taken from Ref. jHj. 
Right: Allowed region at 90%, 95%, 99%, 99.73% (3cr) C.L. for the neutrino oscillation 
parameters Am 2 = Am 2 ol = Am|; amLAND and tan 2 6 = tan 2 i) so \ = tan 2 '^KamLAND ob¬ 
tained in Ref. IBS! from the combined analysis of solar and KamLAND data (see also 

Refs, mmmmmmmm ) • 


liquid scintillator detector of the KamLAND experiment is located in the Kamioka mine, 
in Japan, at a depth of about 1 km. Both prompt photons from the annihilation of e + in 
the scintillator and 2.2 MeV delayed photons from the neutron capture n + p —> d±y are 
detected (the mean neutron capture time is 188 ±23 /rsec). In order to avoid background, 
mainly from the decays of 238 U and 232 Th in the Earth, the cut -E prompt > 2.6 MeV was 
applied. 

During 145.1 days of running 54 events were observed. The number of events expected 
in the case of no neutrino oscillations is 86.8 ± 5.6. The ratio of observed and expected 
u e events is 


Aobs — Ai 


BG 


N ( 


expected 


0.611 ±0.085 ±0.041, 


(2.35) 


where Abg = 0.95 ± 0.99 is the estimated number of background events. 

The left figure in Fig. [3] shows the dependence of the ratio of observed and expected 
v e events on the average distance between reactors and detectors for all reactor neutrino 
experiments. The dotted curve was obtained with the best-fit solar neutrino LMA values 
of the oscillation parameters Am 2 ol = 5.5 x 10~ 5 eV 2 and sin 2 2d so i = 0.83 obtained in 
Ref. IB2J. 

In the KamLAND experiment the prompt energy spectrum was also measured (see 
the right figure in Fig. na>. The prompt energy is connected with the energy of v e by 
the relation E prompt = E Pe — 0.8 MeV — E n ( E n is the average kinetic energy of the 
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neutron and 0.8 MeV = m n — rn p — m e , with the electron neutrino mass coining from the 
annihilation of the final positron in Eq. (EH with an electron in the medium). From 
the two-neutrino analysis of the KamLAND spectrum the following best-fit values of the 
oscillation parameters were obtained: 

^KamLAND = 6 -9 X 10 _5 eV 2 , sill 2 2d K amLAND = 1 ■ (2.36) 

The 95% C.L. allowed regions in the plane of the oscillation parameters obtained from 
the analysis of the measured rate and energy spectrum are shown in the left figure in 
Fig. EH The region outside the solid line is excluded from the rate analysis. The dark 
region is the solar neutrino LMA allowed region, obtained in Ref. One can see that 
two of the KamLAND allowed regions overlap with the solar neutrino LMA region. 

The KamLAND result provides strong evidence of neutrino oscillations, obtained 
for the first time with terrestrial reactor antineutrinos with the initial flux well under 
control. The KamLAND result allows to exclude the SMA, LOW and VAC solutions 
of the solar neutrino problem. It proves that the only viable solution of the problem is 
LMA. The right figure in Fig. El shows the allowed region for the oscillation parameters 
obtained in Ref. J25] from the combined analysis of solar and KamLAND data (see also 

Refs. jHHEHlIHiEniEIliniESlEll)- 

2.4 CHOOZ and Palo Verde 

The results of the long-baseline reactor experiments CHOOZ [M] and Palo Verde [%H] . in 
which D e disappearance due to neutrino oscillations in the atmospheric range of Am 2 was 
searched for, are very important for the issue of neutrino mixing. In these experiments 
electron antineutrinos were detected via the observation of the process 

v e + p —> e + + n . (2.37) 

No indication in favor of the disappearance of reactor P e ’s was found. The ratio R of the 
measured and expected numbers of v e events in the CHOOZ j2| an d the Palo Verde [HSj 
experiments are, respectively, 

R= 1.01 ±2.8% ±2.7%, R= 1.01 ±2.4% ±5.3%. (2.38) 

From the 95% C.L. exclusion plot obtained in Ref. fM] from the two-neutrino analysis 
of CHOOZ data, for Am^nooz = 2-5 x 10 _3 eV 2 (the S-K best-fit value for Am 2 tm , see 
Eq. (12.171) 1 we have 

sin 2 2tf C Hooz < 1-5 x HT 1 . (2.39) 


2.5 Phenomenology 

In the minimal scheme with mixing of three massive neutrinos, the 3x3 Pontecorvo-Maki- 
N akagawa- S akat a 13511331321 mixing matrix U is characterized by three mixing angles 
and one CP phase (in the case of Dirac neutrinos; in the case of Majorana neutrinos, 
in the mixing matrix there are two additional phases which are irrelevant for neutrino 
oscillations). Let us discuss now neutrino oscillations in the atmospheric and solar ranges 
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Figure 7: The two types of three-neutrino mass schemes compatible with the hierarchy 
Am 2 oi < A 


of neutrino mass-squared differences in the framework of this scheme, which provides two 
independent Am 2 ’s: Am^ = m\ — m\ and Am § 2 = m| — m|. 

Two important features of the neutrino mixing, which were revealed in the recent 
solar, atmospheric, long-baseline reactor and accelerator experiments, determine neutrino 
oscillations. 

The first feature is the hierarchy of the neutrino mass-squared differences: from the 
analyses of the data of the solar and atmospheric neutrino experiments it follows that 
Am 2 ol <C Arn 2 tm . This hierarchy can be realized only with the two types of three-neutrino 
mass schemes 7 shown in Fig. [3 The absolute scale of the neutrino masses in the two 
schemes in Fig. |3 is not fixed by neutrino oscillation experiments. Figure |S] shows the 
neutrino masses as functions of the lightest mass rn 1 (see Ref. |9Z])- One can see that 
the “normal” scheme in Fig. |3 is compatible with the natural mass hierarchy (12. 711 if 
mi -C m 2 , whereas in the “inverted” scheme // 2 and 1/3 are always almost degenerate. 

Let us first consider the “normal” mass scheme in Fig. [3 with Am 2 ! — Am 2 ol and 
A)71 j 2 ~ Am 2 tm . In this case we have 

AmJx <C Am .32 . (2.40) 

For the values L/E relevant for neutrino oscillations in the atmospheric range of Am 2 (at¬ 
mospheric and long-baseline accelerator and reactor experiments) we have A m^L/E -C 
1. Hence, we can neglect the contribution of A/n 2 ! to the transition probability in 
Eg. (12.111). For the Drobabilitv of the transition u.. —* u,.y. we obtain (see Ref. ESI) 


P(z/ a 


zv) 



^1 — cos Am. 3 2 



(a a '), 


(2.41) 


7 Independently from the type of neutrino mass spectrum (“normal” or “inverted”), we label neutrino 
masses in such a way that m\ < m 2 < m 3 . Another convention is often used in the literature, such 
that in both “normal” and “inverted” mass spectra Am s 2 Q | ~ A> 0 and Ara 2 tm ~ |Am| 2 |- In this 
notation, mi < m 2 < m 3 in “normal” schemes and m 3 < mi < m 2 in “inverted” schemes. 
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Figure 8 : Neutrino masses as functions of the lightest mass m\ in the two three-neutrino 
mass schemes in Fig. 0 The dashed and dotted lines represent, respectively, the limits 
for m 2 and m 3 from solar and atmospheric neutrino data. 


where the oscillation amplitude is given by 

K'-, a = 4\U a >3\ 2 \U a3 \ 2 (2.42) 


In the standard parametrization of the mixing matrix (see Ref. 


m ) we have 


U e 3 = sin #13 e 16 , = yjl - |£ 43 1 2 sin d 23 , U r3 = \A - \U e3 \ 2 cos #23 , 

(2.43) 

where #13 and #23 are mixing angles and 6 is the CP-violating phase. Hence, for the 
amplitude of the transitions —> v T and v fl —> v e we obtain, respectively, 


A t;m = (1- |f/ e3 | 2 ) 2 sin 2 2#23, A e;/i = 4 \U e3 \ 2 (l — |I/ e3 | 2 ) sin 2 # 23 • 
The probability of v a to survive is given by 


PK 


v n ) = 


1 - E p (^ 

ol'^ol 


V a ') = 1-B a;Q , I 1 — COS 


L \ 


: 


2 E J 


(2.44) 


(2.45) 


where 

B Q-a = = ^ l^asl 2 (l — |Hq, 3 | 2 ) . (2.46) 

Thus, due to the hierarchy in Eq. (12.401) . the transition probabilities in the atmospheric 
range of A?n 2 have a two-neutrino form. Taking into account that the elements |t/ a3 | 2 , 
which determine the oscillation amplitudes, satisfy the unitarity condition I £ 4 * 31 2 = 
1 , we conclude that the transition probabilities are characterized by three parameters: 
A? 77 2 2 , sin 2 2#23 and | C/ e 3 1 2 • 
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Let us consider now neutrino oscillations in the solar range of Am 2 . The v e survival 
probability in vacuum can be written in the form 

2 

J2 | U ei \ 2 e~ lArn ^ + \U e3 \ 2 e~ iAm2 ^ (2.47) 

4 = 1,2 



We are interested in the survival probability averaged over the region where neutrinos 
are produced, over the detector energy resolution, etc.. Because of the neutrino mass- 
squared hierarchy in Eq. flZHD , in the expression for the averaged survival probability 
the interference between the first and second terms in (HUD vanishes and we obtain 


where 


and 


P(U £ - U e ) = |f/ e3 | 4 + (1 - |£4 3 | 2 ) 2 P {1 ’ 2 \v e - »e) , 

p(1.2) 


1 

2 


[v P —> u P ) — 1-A*' 1,2 -’ ( 1 — cos Am7 —r I 


2 L_\ 

21 2 EJ ’ 


(1 - IVal 2 ) 2 

We have used the standard parametrization of the neutrino mixing matrix (see Eq. 
and Ref. JS]), with 


(2.48) 

(2.49) 

(2.50) 

El 


u e 1 = \/l - |t/ e3 | 2 cos di 2 , f/ e2 = \A - IV e3 1 2 sini?i 2 , (2.51) 


where i9i2 is a mixing angle. The expression (I2.48j) is also valid in the case of oscillations 
in matter HH. In this case P( 1,2 )(z/ e —► u e ) is the two-neutrino survival probability in 
matter calculated with the charged-current matter potential Vcc multiplied by (1 —|P e3 | 2 ). 

The second important feature of the neutrino mixing is the smallness of the parameter 
|t/ e3 | J . This follows from the results of the CHOOZ and Palo Verde experiments and from 
the results of solar neutrino experiments. In the CHOOZ and Palo Verde experiments, 
the probability of P e to survive is 

P(p e -> P e ) = 1 - i B e;e fl - cosAm^j , (2.52) 


where 

B e;e = 4|V e3 | 2 (l-|V e3 | 2 ) . (2.53) 

From the two-neutrino exclusion plots, obtained in the CHOOZ and Palo Verde experi¬ 
ments mm, it follows that 

(2.54) 

where the upper bound Bg. e depends on Am 2 2 . 

For the S-K [1,0,3.; allowed values of Am 2 2 , from the 95% C.L. CHOOZ exclusion 
plot we find 

lxlO- 1 <B° ;e <2.4xlO- 1 . (2.55) 
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Thus, the parameter |C/ e 3 1 2 can be small or large (close to one). This last possibility 
is excluded by the solar neutrino data (see DSD- At the S-K best-fit point A m\ 2 = 
2.5 x lCT 3 eV 2 we have 


\u e3 \ 2 < 4 x 10" 2 

A combined fit of all data leads to [BBJ 

(95% C.L.). 

(2.56) 

\U e3 \ 2 < 5 x IQ’ 2 

(99.73% C.L.). 

(2.57) 


There are three major consequences of the neutrino mass-squared hierarchy mi) 
and of the smallness of | t/ e 3 1 2 : 

1. The dominant transition in the atmospheric range of Am 2 is —> v T . From 

Eq. (12.4411 it follows that 

Am 2 32 ~ Am 2 tm , sin 2 2^ 23 ~ sin 2 2-d atm . (2.58) 

2. In the solar range of Am 2 the probability of v e to survive has the two-neutrino form 

P(v e ^v e )~pW(v e ^v e ). (2.59) 

From Eqs. m and it follows that 

A/77,21 — Am 2 ol , tan 2 $ 12 ~ tan 2 'd sol . (2.60) 

3. Neutrino oscillations in the atmospheric and solar ranges of Am 2 in the leading 

approximation are decoupled mm- Oscillations in both ranges are described by 
two-neutrino formulas, which are characterized, respectively, by the parameters 
Am| 2 , sin 2 2 d 23 and Am 21 , tan 2 di 2 . 

We have considered the hierarchy of the neutrino mass-squared differences, 

which is realized in the “normal” mass scheme in Fig. 0 The data of neutrino oscillation 
experiments are compatible also with the “inverted” mass scheme in Fig.0 with A m 2l ~ 
Am 2 tm and Am| 2 — Am 2 ol . In this case an inverted hierarchy of the mass-squared 
differences takes place 8 : 

Am\ 2 <C A717-23 (2.61) 

The expressions for the transition probabilities in the case of the inverted hierarchy 
can be obtained from Eqs. (12.411) . (12.421) . (12.451) . (|2.46D . (I2.48I) - (I2.5()I) with the change 
Am| 2 ^ Am 2 !, \U a3 \ 2 \U al \ 2 , 0 12 t? 23 , P™ - P^\ A (1,2) A (2 ’ 3) . 

We have discussed up to now evidences in favor of neutrino oscillations that have been 
obtained in the atmospheric and solar neutrino experiments. There exist at present also 
an indication in favor of short-baseline —> v e transitions, which has been obtained only 
in the accelerator experiment LSND 133- The LSND data can be explained by neutrino 

8 The type of the neutrino mass spectrum (“normal” or “inverted”) may be determined via the inves¬ 
tigation of v e —> Un and v e —> oscillations in future long-baseline experiments if | C/ e 3 1 2 is not too small 
(see Ref. E3). The distance between the neutrino source and the detector in such experiments must be 
large enough for the matter effects to be sizable. 
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oscillations. From a two-neutrino analysis of the data, the best-fit values of the oscillation 
parameters are 

Atttlsnd — 1-2 eV 2 , sin 2 2d LSND — 3 x 1CT 3 . (2.62) 

In order to describe the results of the solar, atmospheric and LSND experiments, 
which require three different values of neutrino mass-squared differences Am 2 ol , A/n^ tm 
and Amj sND, it is necessary to assume that at least one sterile neutrino exists in addition 
to the three active neutrinos z/ e , z/^, v T . In the mass basis, in addition to the three 
light neutrinos V \, z/ 2 , rq there must be at least one neutrino with mass of the order 
a/A mL SND ~ 1 eV (see, for example, Ref. [TTJ]). However, in spite of the additional degrees 
of freedom, schemes with four neutrinos do not fit well the data (see Refs. (io2inns!) 9 . 

The result of the LSND experiment requires, however, confirmation. The MiniBooNE 
experiment at Fcrmilab f'(T4l . that started recently, is aimed to check the LSND result. 

From neutrino oscillation experiments we can obtain information only on the neutrino 
mass-squared differences, not on the absolute values of neutrino masses. The great advan¬ 
tage of neutrino oscillations experiments, that was stressed in the early papers on neutrino 
oscillations is that they are sensitive to very small values of Am 2 . This is 

connected with the fact that neutrino oscillations are an interference phenomenon. It is 
also important that there is the possibility to perform experiments with detectors at very 
large distances from neutrino sources (solar, atmospheric and long-baseline experiments) 
and for small neutrino energies (solar and reactor experiments). 

The understanding of the origin of neutrino masses and neutrino mixing requires 
knowledge of the absolute values of neutrino masses (see Refs. mm)- The problem 
of the absolute values of neutrino masses is one of the most challenging problems of the 
physics of massive and mixed neutrinos. At present there are only upper bounds for the 
absolute values of neutrino masses. The most stringent bound was obtained from the 
experiments on the measurement of the high energy part of the /3-spectrum of 3 H. In the 
next section we will discuss the results of these experiments and future prospects. 

3 Neutrino mass from ( 3 -decay experiments 

3.1 The /3-spectrum in the case of neutrino mixing 

The method of measurement of the neutrino mass through the detailed investigation of 
the high-energy part of the /3-spectrum was proposed in 1934 by Fermi in his classical 
paper on the theory of /3-decay (1061 and by Perrin BUZ!- The first experiments on the 
measurement of the neutrino mass with this method have been done in 1948-49 mm- 

Usually, the neutrino mass is measured through the measurement of the high energy 
part of the /3-spectrum of tritium 

3 H —> 3 He + e~ + u e . (3.1) 

The investigation of this decay has several advantages. Since tritium decay is superal- 
lowed, the nuclear matrix element is a constant and the electron spectrum is determined 

9 Since there is no experimental indication in favor of transitions of active neutrinos into sterile states, 
schemes with more than four neutrinos are disfavored as well. 
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by the phase space. It has a relatively small energy release Eq ~ 18.6 keV and a conve¬ 
nient lifetime (Tf/ 2 — 12.3 years). A small value of Eq is convenient because the relative 
fraction of events in the high energy part of the spectrum, which is sensitive to the 
neutrino mass, is proportional to Eq 3 . 

Let us consider the decay EH in the case of nonzero neutrino masses and neutrino 
mixing. The effective Hamiltonian of the process is 

7-(-? C = 2 e L ^ a u eL j a + h.c., (3.2) 

where j a is the hadronic charged current and the field v e L is given by (see Eq. (I2.2J) ) 

^eL ^ ^ . (3.3) 

i 


The state of the final particles in the decay EH) is 

I/) = TI P') |p> I Pi) Ki (P,P„ P'KS - l)l^>. (3.4) 

i 

Here p is the momentum of the electron, P and P' are the momenta of the initial and 
final nuclei, pi is the antineutrino momentum, | P'), | p), \pi) are the normalized states of 
the final nucleus, electron, antineutrino with mass mj, and 

(p,Pi,P'\(S - 1)1 P) = (p,Pi, P'\T\P) (2vr f5\ P + Pi + P’-P) (3.5) 


is the element of the A-matrix. 

We are interested in the spectrum of electrons. After the integration over the momenta 
P', pi and over the angle of emission of the electron, we have 



i 


(3.6) 


where 


-^ = Cp(E + m e ) (E 0 - E) \J(E 0 — E) 2 — F(E) 6(E 0 - E - m % ) . (3.7) 

Here E is the kinetic energy of the electron, Eq is the energy released in the decay, m e 
is the mass of the electron and F{E) is the Fermi function, which describes the Coulomb 
interaction of the final particles. The constant C is given by 

e; 

Ttl 

C = G 2 f —-^ cos 2 0 c \M\ 2 , (3.8) 

2 7T- 3 

where Gf is the Fermi constant, Qq is the Cabibbo angle, M is the nuclear matrix element 
(a constant). 

Let us notice that neutrino masses enter in the expression EH through the neutrino 
momentum \p t \ = \J ( Eq — E) 2 — rrif. The step function 6(E 0 — E — rrii ) provides the 
condition E < Eq — rrii. The recoil of the final nucleus was neglected in Eq. EH- 
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Two experiments on the measurement of neutrino masses with the tritium method 
are going on at present (Mainz [3T11321133] and Troitsk [331133] ). The sensitivity of these 
experiments to the neutrino mass is about 2-3 eV. The sensitivity to the neutrino mass of 
the future experiment KATRIN |3Sj is expected to be about one order of magnitude better 
(0.35 eV). We will discuss the results of the Mainz and Troitsk experiments later. Now 
we consider the possibility to determine the neutrino mass from the results of the /3-decay 
experiments for different neutrino mass spectra, having in mind these sensitivities. 

As it is seen from Eq. (EH). the largest distortion of the /3-spectrum due to neutrino 
masses can be observed in the region 


E 0 — E ~ rrti. (3.9) 

However, for m l ~ leV only a very small part (about 1CU 13 ) of the decays of tritium 
give contribution to the region El. This is the reason why in the analysis of the results 
of the measurement of the tritium /3-spectrum a relatively large part of the spectrum is 
used (for example, in the Mainz experiment the last 70 eV of the spectrum). Taking this 
into account, the tritium /3-spectrum that is used for the fit of the data can be presented 
in the form (see also the discussion in Ref. [114] ) 

= Cp {E + m e ) {Eq - E) ^J{EQ-Ey-m}F{E) , (3.10) 

where the effective mass rrig is given by 

m 0 = Yl\ u <*\ 2m ? ( 3 - n ) 

i 

Let us consider first the minimal scheme with three massive neutrinos iq, zq and v 3 . 
The minimal neutrino mass m\ and the character (“normal” or “inverted”, hierarchical 
or almost degenerate) of the neutrino mass spectrum are unknown at present. Neutrino 
oscillation experiments allow to determine the neutrino mass-squared differences 
and Am| 2 - Hence, it is possible to express the values of the neutrino masses m 2 and m 3 
in terms of the unknown mass 777 1 as 


771 2 = \]m\ + A m%i , m 3 = \Jm\ + Am^ + A m \ 2 . (3-12) 

In the “normal” three-neutrino scheme in Fig. [TJ Am^ = Am( ol and A 77732 = A/77^ tm . 
Using Eqs. (12311) and (123(1 . we obtain 

m 2 g = ml + (sin 2 t? sol + cos 2 t9 so i |U e3 | 2 ) Am 2 ol + |U e3 | 2 Am 2 tm . (3.13) 

I 11 the case of the natural neutrino mass hierarchy (EH), we have 

7772 — \J A/77 g ol ~ 7 X 10 -3 eV , 7773 ~ \J ATTl^tm — 5 X 10 _2 eV , (3.14) 

where the best-fit values (TT71) and (13331) of the oscillation parameters were used. 

The contribution of the heaviest neutrino mass m 3 to the effective neutrino mass 
Emu enters with the weight |I/ e3 | 2 , for which we have the upper bound (12.571) from the 
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results of the CHOOZ experiment JEH]. Taking into account this bound and using the 
best-fit value in Eq. dznnp for 'dsoi, for the effective neutrino mass m@ we obtain 

mp ~ (sin 2 $ sol Am 2 ol + \U e3 \ 2 A/n 2 tm ) 1/2 < 1.2 x 10~ 2 eV , (3.15) 

which is about one order of magnitude smaller than the sensitivity of the future tritium 
experiment KATRIN j3S]. 

In the “inverted” three-neutrino scheme in Fig. H with A= Am 2 tm and Am | 2 — 
Am 2 ol , using Eq. (12.5111 . we obtain 

m 2 g = ml+(l- \U el \ 2 ) (Arn 2 tm + sin 2 d sol Am 2 ol ) . (3.16) 

The value of \U e \\ 2 in the “inverted” neutrino scheme is bounded by the results of the 
CHOOZ experiment as the value of \U e3 \ 2 in the “normal” neutrino scheme (see Eq. (12.5711 
and the remark after Eq. (HMD): 

\U e il 2 < 5 x 10 ” 2 (99.73% C.L.). (3.17) 

In the case of the “inverted” neutrino mass hierarchy m\ -C m 2 < m 3 we have 

m 2 ~ m 3 ~ \/ A m 2 tm ~5x 1CT 2 eV , mi< a/ Am 2 tm . (3.18) 

For the effective neutrino mass mp we obtain 

mp ~ a/A m 2 tm — 5 x 1CT 2 eV , (3.19) 

which is also much smaller than the sensitivity of the KATRIN experiment. 

Figure El shows the allowed values of mp as a function of m\. One can see that in both 
the “normal” and “inverted” three-neutrino schemes in Fig. [TJ the KATRIN experiment 
may obtain a positive result only if the three neutrino masses are almost degenerate and 
mi is of the same order or larger than the sensitivity of the experiment (mi > 0.3eV). 
In this case m\ ~ m 2 — m 3 , and from the unitarity of the mixing matrix we obtain 

mp ~ mi, (3.20) 


as shown in Fig. El 

If the LSND result [21; is confirmed by the MiniBooNE [ TIM; experiment, it will mean 
that (at least) four massive and mixed neutrinos exist in nature 10 . 

Let us discuss now the possibilities to measure the neutrino mass with the tritium 
method in the case of four massive neutrinos. In this case, we have three different neutrino 
mass-squared differences Am 2 ol , Am 2 tm and Am 2 SND , given by (12.1711 . (|2.33j) and (12.6211 . 
Let us assume that mi -C \/ Am^ND- 

Figure HTH shows the six four-neutrino mass spectra compatible with the mass-squared 
hierarchy Am 2 ol -C Am 2 tm -C A?71 lsnd- I n all spectra there are two groups of close 
masses, separated by the LSND gap of the order of 1 eV. There are two possibilities for 
the groups: 2+2 and 3+1. 

10 Although the relatively bad fit of the data in the framework of four-neutrino schemes (see Refs. m 
GUI) may suggest the possibility of more exotic explanations. 
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Figure 9: Value of mg as a function of rri j (see Eqs. (13.FID and (13.1 (ID ). The two dashed 
curves have been calculated assuming the best-fit values of Am 2 ol and i ? so i in Eq. (El 
and the best-fit value of Am 2 tm in Eq. (hud. The lower dashed curves correspond 
to |C/ e 3 1" = 0 in the “normal” scheme and \U e i\ 2 = 0 in the “inverted” scheme. The 
upper dashed curves correspond to the upper limits | C/ e 3 1 2 = 5 x ICE 2 (Eq. (12.57(1 ) in the 
“normal” scheme and \U e i\ 2 = 5 x ICE 2 (Eq. (13 .1 711 1 in the “normal” scheme. The two 
solid curves represent the lower and upper limits for mg obtained from the 99.73% C.L. 
solar LMA-MSW region in Ref. and the 99% C.L. atmospheric region in Ref. da, 
and \U e3 \ 2 and |f/ e i| 2 bounded by Eqs. (I2.57D and (I3.17D . respectively, in the “normal” 
and “inverted” schemes. 

In order to calculate the contribution of neutrino masses to the /9-spectrum it is 
necessary to take into account the constraints on the elements of the neutrino mixing 
matrix that can be obtained from the data of the short-baseline reactor experiment Bugey 
[11(3] . in which no indication in favor of neutrino oscillations was found. 

In the framework of four-neutrino mixing, the probability of the reactor h e ’s to survive 
is given by (see fT^j l 


P^e 




^1 — cosAm^ 



(3.21) 


Here 



where i runs over the mass indices of the first (or second) group. 

From the exclusion curve obtained in the Bugey experiment [ITBj . we have 


(3.22) 


B P . P < B 


— e;e> 


where the upper bound B° e depends on A 


(3.23) 
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Figure 10: The six allowed types of four-neutrino mass schemes. 


Let us consider first the neutrino mass spectra I, II and B in Fig. lTUl in which the solar 
neutrino mass-squared difference Am;? ol belongs to the light group. Taking into account 
the solar neutrino data, from the Bugey exclusion plot we obtain 

]T|f/ eJ | 2 <2xlO- 2 , (3.24) 

i 

where i runs over indices of neutrinos belonging to the heavy group (i — 3,4 for 2+2 
scheme and i — 4 for 3+1 schemes). Thus, in this case the contribution of heavy neutrinos 
to the (3 spectrum is suppressed. For the effective neutrino mass we have the upper bound 


< 



(3.25) 


Taking into account ()3.24D and using the best fit values of the parameters AmL SND and 
Am^tm (see (12.171) and (12.(121) ). for the effective neutrino mass mp we obtain the bound 


mp < 1.6 x 10 1 eV, 


(3.26) 


which is smaller than the sensitivity of the future tritium experiment KATRIN. 

If Ar7ig 0l is the mass-squared difference of neutrinos belonging to the heavy group, in 
Eq. (13.221) the index i takes the values i — 1,2 for 2+2 scheme and i — 1 for 3+1 schemes. 
In this case, taking into account the unitarity of the mixing matrix, we have 

nip ~ yj Amj( SND . (3.27) 

The allowed range for the parameter Am^ND found in Ref. EU is (see also Ref. HE) 

0-2 eV 2 < Am£ SND < 2 eV 2 . (3.28) 

From Eqs. H23) and (13.281) . for the effective neutrino mass we have the limits 


0.45 eV < nig < 1.4 eV . 


(3.29) 


In this case, if the data of the LSND experiment are confirmed by the MiniBooNE 
experiment, there is a chance to observe the effect of neutrino mass in the future KATRIN 
experiment [3S with the expected sensitivity of about 0.35 eV. 
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Figure 11: Left: Integral spectrum measured in 1994 (open circles), 1998/1999 (filled 
squares) with fit (line) and 2001 (open squares) in the Mainz experiment. Right: Mainz 
fit results on ml = as a function of the lower boundary E\ ow of the fit interval for 
data from 1998 and 1999 (open circles) and from 2001 (filled circles). Figures taken from 
Ref. [3H] , 


3.2 Mainz experiment 

The source in the Mainz experiment is frozen molecular tritium condensed on 

a graphite substrate. The spectrum of the electron is measured by an integral MAC-E- 
Filter spectrometer (Magnetic Adiabatic Collimator with a retarding Electrostatic filter). 
This spectrometer combines high luminosity with high resolution. The resolution of the 
spectrometer is 4.8 eV. In the analysis of the experimental data four variable parameters 
are used: the normalization C , the background B, the released energy E$ and the effective 
neutrino mass-squared m 2 p. From the fit of the data it was found that E 0 = 18.575 keV. 

The left figure in Fig. [TT] shows the integral spectrum, measured in 1994, in 1998- 
1999, and in 2001, as a function of the retarding energy near the endpoint E 0l and the 
effective endpoint Eo >e ff. The position of E 0)B s takes into account the width of the response 
function of the setup and the mean rotation-vibration excitation energy of the electronic 
ground state of the ^HeT" 1 " daughter molecule. The solid curve was obtained from the fit 
of the data under the assumption mp = 0. Different parts of the spectrum were used in 
the analysis of the data. The right figure in Fig. HT1 shows the dependence of on the 
lower limit Ei ow of the corresponding fit interval (the upper limit is fixed at 18.66 keV, 
well above Eq) for data from 1998 and 1999 (open circles) and from the last runs of 2001 
(filled circles). The error bars show the statistical uncertainties (inner bar) and the total 
uncertainty (outer bar). The correlation of data points for large fit intervals is due to 
the uncertainties of the systematic corrections, which are dominant for fit intervals with 
a lower boundary E\ ow < 18.5 keV. 

In the last 70 eV of the spectrum the combined statistical and systematical error is 
minimal. From the fit of the 1998-1999 experimental data in this interval, it was found 
that 

rrVp = (—1.6 ± 2.5 ± 2.1) eV 2 , (3.30) 


27 





















which corresponds to the upper bound 


171/3 < 2.2 eV (95% C.L.). (3.31) 

In 2001 additional measurements of the /Tdecay spectrum were carried on by the Mainz 
group. From the analysis of these data it was obtained [33] 

m 2 g = (0.1 ± 4.2 ± 2.0) eV 2 . (3.32) 

From the combined analysis of 1998, 1999 and 2001 data it was found [331 

m 2 g = (—1.2 ± 2.2 ± 2.1) eV 2 . (3.33) 

This value corresponds again to the upper bound [231 

mg < 2.2 eV (95% C.L.), (3.34) 

showing that the Mainz experiment has reached his sensitivity limit. 

3.3 Troitsk experiment 

In the Troitsk neutrino experiment [21ll23j, as in the Mainz experiment, an integral elec¬ 
trostatic spectrometer with a strong inhomogeneous magnetic field, focusing the electrons, 
is used. The resolution of the spectrometer is 3.5 — 4 eV. An important difference between 
the Troitsk and Mainz experiments is that in the Troitsk experiment the tritium source is 
a gaseous molecular source. Such a source has important advantages in comparison with 
the frozen source: there is no backscattering, there are no effects of the self-charging, 
the interaction between tritium molecules can be neglected, etc.. In the analysis of the 
data the same four variable parameters C, B , E 0 and m 2 g were used. From the fit of 
the data, for the parameter m 2 g large negative values in the range —20 -j—10 eV 2 have 
been obtained. The investigation of the character of the measured spectrum suggests 
that the negative m 2 g is due to a step function superimposed on the integral continuous 
spectrum. The step function in the integral spectrum corresponds to a narrow peak in 
the differential spectrum. 

In order to describe the data, the authors of the Troitsk experiment added to the 
theoretical integral spectrum a step function with two additional variable parameters 
(position of the step E step and the height of the step). From a six-parameter fit of the 
data, the Troitsk Collaboration found 

m 2 g = (—2.3 ± 2.5 ± 2.0) eV 2 . (3.35) 

This value corresponds to the upper bound 

mg < 2.2 eV (95% C.L.), (3.36) 

as in the Mainz experiment. 

The Troitsk Collaboration found that the position of the step E 0 — E step changes 
periodically in the interval 5 — 15 eV and that the average value of the height of the step 
is about 6 x ICC 11 of the total number of events. This effect has been called “Troitsk 
anomaly”. Since the Mainz data do not show any indication of a Troitsk-like anomaly, it 
is believed [22 that the Troitsk anomaly is caused by some experimental artifact. 


3.4 Other experiments 

We have discussed up to now tritium experiments for the measurement of the neutrino 
mass. The groups in Genova I1T81 and Milano ma are developing low temperature 
cryogenic detectors for the measurement of the /3-decay spectrum of 187 Re. This element 
has the lowest known energy release (E 0 = 2.5 keV). The relative fraction of events in 
the high energy part of the spectrum is proportional to Eq 3 . Thus, decays with low 
Eq values are very suitable for calorimetric experiments in which the full spectrum is 
measured. The limit for the neutrino mass that was obtained by the Genova group in 
Ref. HHI is 

mp < 26 eV (95%C.L.). (3.37) 

In the future a sensitivity of about 10 eV is expected to be reached. 

3.5 The future KATRIN experiment 

In the future tritium experiment KATRIN ra, two tritium sources will be used: a gaseous 
molecular source (T 2 ), as in the Troitsk experiment, and a frozen tritium source, as in the 
Mainz experiment. The windowless gaseous tritium source will allow to reach a column 
density 5 x 10 1 ' T 2 /cm 2 . 

The integral MAC-E-Filter spectrometer will have two parts: the pre-spectrometer, 
which will select electrons in the last part (about 100 eV) of the spectrum, and the main 
spectrometer. This spectrometer will have a resolution of 1 eV. 

It is planned that the KATRIN experiment will start to collect data in 2007. After 
three years of running the accuracy in the measurement of the parameter m 2 g will be 
0.08eV 2 . This will allow to reach a sensitivity of 0.35 eV in the determination of the 
effective neutrino mass mp. 

As in the case of the Mainz experiment, in the analysis of the data of the future 
KATRIN experiment four variable parameters are planned to be used. The value of the 
parameter E 0 can be taken, however, from the independent measurement of the 3 H and 
3 He mass difference. If the accuracy of such measurements reaches 1 p.p.m., the sensitivity 
of the KATRIN experiment to the neutrino mass will be signihcantly improved. 

In the KATRIN experiment, not only the integral spectrum, but also the differen¬ 
tial spectrum is planned to be measured. These measurements will allow to clarify the 
problem of the Troitsk anomaly in a direct way. 

4 Muon and tau neutrino mass measurements 

Information on the “mass” of the muon neutrino can be obtained from the measurement 
of the muon momentum in the decay 

vr + —> n + + Uf ,. (4.1) 

We discuss here such measurements from the point of view of neutrino mixing. From 
four-momentum conservation in the decay HUD it follows that 

= m l + m l~ 2 "G yjml + p 2 . (4.2) 
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Here m, is the mass of z/,, and m M are the masses of the pion and the muon, and p M 
is the muon momentum in the pion rest-frame. 

The value of the muon momentum measured in the most precise PSI experiment Q23! 

is 

p= 29.79200 ± 0.00011 MeV. (4.3) 

Taking into account the resolution in the measurement of the momentum of the muon 
and the values of the neutrino mass-squared differences measured in neutrino oscillation 
experiments (see Section EJ), we come to the conclusion that the effect of neutrino masses 
in the decay EH can be observed only in the case of an almost degenerate neutrino mass 
spectrum with m 1 3> \J A ml tm (or mi 3> \J A?7Tl SND in the case of four neutrinos). 

In this case, from the unitarity condition JA U IM \ 2 = 1 it follows that the experiments 
on the measurement of the momentum of the muon produced in the decay EH) allow to 
obtain information about the mass rri\ . The value of mi found in Ref. H23! is 

m\ = -0.016 ± 0.023 MeV 2 (vr + -> p+z^). (4.4) 

For the masses of the muon and pion the following values were used » 

= 105.658389 ± 0.000034 MeV, (4.5) 

m n = 139.56995 ± 0.00037 MeV. (4.6) 

The upper bound of the neutrino mass given by the Particle Data Group [031 i s 

mi < 190keV (90%C.L.) (tt+ -> fjt + v„). (4.7) 

The most stringent upper bound on the “mass” of the tau neutrino was obtained in 
the ALEPH experiment 121 .. In this experiment the decays 

T~ —> 27T - 7T + V t , —> 37T” 27r + (7T°) v T (4-8) 

were studied. From the conservation of the four-momentum in the decay 


nir + is T 


(4.9) 


we have 


K = 


m 2 + ml 


rm 


2 m, 


(4.10) 


Here m T is the tau mass, is the invariant mass of the n pions and is the total 
energy of the n pions in the rest frame of the tau. Information on the neutrino mass was 
obtained in Ref. czn from the fit of the distribution d 2 T/dE^dm^ (E h is the total energy 
of the pions in the laboratory system). In the case of neutrino mixing, information about 
the minimal (common) mass mj of an almost degenerate neutrino mass spectrum can be 
obtained from such experiments. The bound obtained in the ALEPH experiment HZP is 


mi < 18.2 MeV (95%C.L.) (r —► mr + u T ). (4-11) 

Thus, the experiments on the measurement of the muon momentum in the decays 
of pions and the d 2 T/dEhdmh distribution in the decays of taus are much less sensitive 
to the absolute neutrino mass mi than tritium experiments. These experiments could, 
however, reveal the existence of particles with masses much larger than the light neutrino 
masses. 


30 





5 Neutrinoless double-/? decay 

The search for neutrinoless double-/? decay 

(A, Z ) —■» ( A , Z + 2) + e + e (5-1) 


of some even-even nuclei is the most sensitive and direct way of investigation of the nature 
of the neutrinos with definite masses (Majorana or Dirac?). Neutrinoless double-/? decay 
is allowed only if massive neutrinos v t are Majorana particles. 

We will assume that the Hamiltonian of the process has the standard form in Eq. (E2) 
and the flavor field v e p is given by the relation (see Eq. (E2 1) 

^eL ^ ^ U ei l^iL t ( 5 . 2 ) 

i 

where i\ are Majorana fields which satisfy the condition 

Vi = vl = C vf . (5.3) 

Here C is the charge conjugation matrix (Cq^tD 1 = —q a , C T = —C). 

The neutrinoless double-/? decay ((/?/?)o^-decay) is a process of second order in the 
Fermi constant Gp, with virtual neutrinos. In the case of the Majorana neutrino mixing 
in Eq. Q, the neutrino propagator is given by the expression 

( 0 | T(v eL ( Xl )^ L (x 2 m ^ (m) I ZpLc. ( 5 . 4 ) 

Here 

W = U ei m i ■ ( 5 -5) 

i 

The matrix element of neutrinoless double-/? decay is proportional to the nuclear matrix 
element and to the effective Majorana mass (m), which depends on neutrino masses m* 
and on U^. 

The elements of the neutrino mixing matrix U e i are complex quantities. In the case 
of CP invariance in the lepton sector, the elements U e i satisfy the condition [mrm 

u* ei = v* Uei , (5.6) 


where rji = i pi is the CP-parity of the neutrino i\ (/q = ±1). Let us write down 


U ei 


\U e , 




(5.7) 


From Eq. (E3 we obtain 

= (5-8) 

Thus, in the case of CP invariance in the lepton sector, the effective Majorana mass is 
given by 

(m) = ^2\Uei \ 2 e l % Pi rrii . (5.9) 
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Experiment 

Nucleus 

Sensitivity 

h°r 2 [yr]' 

Sensitivity 

IWI [eV] 

NEMO 3 P3Dj 

100 Mo 

4 x 10 24 

5.6 x 10" 1 

COBRA PSU 

130 Te 

1 x 10 24 

2.4 x 10" 1 

CUORICINO H323 

130 Te 

1.5 x 10 25 

1.9 x 10" 1 

XMASS H33j 

136 Xe 

3.3 x 10 26 

9 x 10~ 2 

CAMEO mm 

116 Cd 

1 x 10 26 

6.9 x 1(T 2 

exo m 

136 Xe 

8 x 10 26 

5.2 x 1(T 2 

MOON TiMCTT?) 

100 Mo 

1 x 10 27 

3.6 x HT 2 

CUORE H32J 

130 Te 

7 x 10 26 

2.7 x 1(T 2 

Majorana [iTSHlj 

76 Ge 

4 x 10 27 

2.5 x HT 2 

GEM mi 

76 Ge 

7 x 10 27 

1.8 x 1(T 2 

GENIUS HD] 

76 Ge 

1 x 10 28 

1.5 x HT 2 


Table 1: Future neutrinoless double-/? decay projects. 

The results of many experiments on the search for (P/3) o^-decay are available at present 
(see Refs. mm)- No indication in favor of (/?/3) 0l/ -decay have been obtained up to 
now 11 . The most stringent lower bounds for the lifetime of (/?/?) 0l ,-decay were obtained 


in the Heidelberg-Moscow [12.7J and IGEX | l28j 76 Ge experiments: 

Ty 2 > 1-9 x 10 25 yr (90%C.L.) (Heidelberg-Moscow), (5.10) 

T,°; 2 > 1.57 x 10 25 vr (90%C.L.) (IGEX). (5.11) 

Taking into account different calculations of the nuclear matrix element, from these results 
the following upper bounds for the effective Majorana mass were obtained: 

\(m)\ < (0.35 — 1.24) eV (Heidelberg-Moscow), (5.12) 

|(m)| < (0.33- 1.35) eV (IGEX). (5.13) 

Many new experiments for the search of neutrinoless double-/? decay are in preparation 
at present (see Tableland Ref. j22])- In these experiments the sensitivities 

| (m) | ~ 1.5 x HT 2 - 5.6 x HT 1 eV (5.14) 


are expected to be achieved 12 . 

An evidence for neutrinoless double-/? decay would be a proof that neutrinos with 
definite masses are Majorana particles and that neutrino masses have an origin beyond 
the Standard Model. The value of the effective Majorana mass \ (m) | combined with the 
results of neutrino oscillation experiments could allow to obtain important information 
about the character of the neutrino mass spectrum, about the minimal neutrino mass mi 
and about the Majorana CP phase (see Refs. prrirraflrrasin 44in 45irran] and references 
therein). 

11 The recent claim ED of an evidence of the (/?/?) o^-decay, obtained from the reanalysis of the data 
of the Heidelberg-Moscow experiment, was strongly criticized in Refs. rmrm 

12 These sensitivities have been estimated using the nuclear matrix elements calculated in Ref. cm 
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Let us consider three typical neutrino mass spectra in the case of three massive and 
mixed neutrinos 13 : 


1 . m i m 2 m 3 (hierarchy of neutrino masses). 

For the effective Majorana mass \(m)\ we have the upper bound 

\(m) \ < sin 2 $ sol Am 2 sol +\U e3 \ 2 x/AmLn- (5-15) 

Using the best-fit values of the oscillation parameters Am 2 tm , Am( ol , tan 2 , d so i, and 
the CHOOZ bound on \U e3 \ 2 (Eqs. (12.171) . (12.331) . and (12.561) 1. we have 

|(m)| < 3.8 x 10“ 3 eV. (5.16) 


Taking into account the upper bounds for the oscillation parameters, one obtains 
|148j 

\{m)\ < 8.2 x 10- 3 eV. (5.17) 

The bounds (15.161) and flUZD are signihcantly smaller than the expected sensitivities 
(EH of the future (/3/3) 0u -decsLy experiments. Thus, the observation of ((3/3) 0u - 
decay in the experiments of the next generation could pose a problem for the natural 
hierarchy of neutrino masses. 

2. mi -C m 2 < m% (inverted hierarchy of neutrino masses). 

The effective Majorana mass is given by 

|(m)| ~ (l - sin 2 2't? sol sin 2 a ) 5 a/A m 2 tm , (5.18) 


where a = 0:3 — «2 is the the difference of CP phases. From this expression it follows 
that 

V Am 2 tm I cos 2tf sol | < I (m) | < \/ Am 2 tm , (5.19) 

where the upper and lower bounds correspond to equal and opposite CP parities in 
the case of CP conservation. 

Using the best-fit value of the parameter tan 2 in Eq. (EH) , we have 

\ V 7 Am 2 tm < | (m) | < v 7 Am 2 tm . (5.20) 


Thus, in the case of an inverted mass hierarchy, the scale of |(m)| is determined 
by a/A m 2 tm ~ 5 x 10~ 2 eV. If the value of \(m)\ is in the range (15.201) . which can 
be reached in the future experiments searching for (/3/3) 0u - decay, it would be an 
argument in favor of an inverted neutrino mass hierarchy. 

The measurement of the effective Majorana mass |(m)| could allow to obtain infor¬ 
mation about the CP phase a mm- Indeed, from Eq. (EmHl we have 


sin' a ~ ( 1 


IWP ( 

sin 2 2 'd so i 


(5.21) 


13 Neutrinoless double-/? decay in the case of four-neutrino mixing was considered in detail in Ref. imi- 
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3. mi — m 2 — m% (almost degenerate neutrino masses). 


Let us assume that mi In this case m 2 — m 3 ~ mi in both the 

“normal” and “inverted” spectra in Fig. Q For the effective Majorana mass, inde¬ 
pendently on the character of the mass spectrum, we have 


| (m) | ~ m-j 



(5.22) 


Neglecting small contribution of |C/ e 3 1 2 (| t/ e 1 ] 2 in the case of the inverted hierarchy), 
for | (m) | we obtain the relations JEEM52U, in which a/A m 2 tm must be changed 
by mi. Thus, if it happens that \(m)\ \J A?n 2 tm ~ 5 x 10 -2 eV, it would be 

a signature of an almost degenerate neutrino mass spectrum. In this case, the 
neutrino mass mi is limited by 

|(m)|<m i < T -L^L< 2 |(m)|. (5.23) 

bl)b Zj U gol 

For the parameter sin 2 a, which characterizes the violation of CP invariance in the 
lepton sector, we have rragirras 1 . 


sm a ~ 



1 

sill 2 2 $ so i 


(5.24) 


If the mass mi is measured in the KATRIN experiment and the precise value 
of the parameter sin 2 219301 is determined in the solar, KamLAND f|, BOREXINO 
G33 and other neutrino experiments, information on the Majorana CP phase can 
be inferred from the results of the future (/9/?)o^-decay experiments. 

Figure H21 [HBJ shows the dependence of |(m)| on mi in the case of the LMA-MSW 
solution of the solar neutrino problem (99.73% C.L. region in Ref. |5[), for the “normal” 
scheme in Fig. Q (left panels) and for the “inverted” scheme in Fig. [3 (right panels). For 
the “normal” scheme with Am 2 ol = Am^, in the case of CP-conservation the allowed 
values of |(m)| are constrained to lie in the medium-gray regions a) between the two 
thick solid lines if 1/21 = ^31 = 1 , b) between the two long-dashed lines and the axes 
if TJ 21 = —?/ 3 i = 1 , c) between the dash-dotted lines and the axes if 1/21 = —7731 = — 1 , 
d) between the short-dashed lines if r/ 2 i = 1/31 = —1. For the “inverted” scheme with 
Am 2 ol = A/TJ 32 , in the case of CP-conservation the allowed regions for |(m)| correspond: 
for \U e i \ 2 = 0.005 and \U e i \ 2 = 0.01 to the medium-gray regions a) between the solid lines 
if rj 2 i = 7731 = ±1, b) between the dashed lines if 1721 = —1731 = ±1; for \U e i \ 2 = 0.05 
to the medium-gray regions c) between the solid lines if 7721 = 7731 = 1 , d) between the 
long-dashed lines if rj 2 i = 7731 = — 1 , e) between the dashed-dotted lines if r] 2 i = —7731 = 1 , 
f) between the short-dashed lines if 7721 = —7731 = — 1. Here 77 ^ is the relative CP-parity 
the neutrinos v l and Uj, given by 

Vij = e ^(pi-pj ). (5.25) 


In the case of CP-violation, the allowed area for |(m)| covers all the gray regions in 
Fig.rnzi Values of | (m) \ in the dark gray regions signal CP-violation. 
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Figure 12: The dependence of |(m)| on mi in the case of the LMA-MSW solution of 
the solar neutrino problem |8j (99.73% C.L.), for the “normal” scheme in Fig. 0 (left 
panels), with Am 2 ol = A/n^i, and for the “inverted” scheme in Fig. 0 (right panels), with 
Am 2 ol = A/n| 2 . Figure taken from Ref. JUD¬ 
AH previous conclusions are based on the assumption that the value of the effective 
Majorana mass |(m)| can be obtained from the measurement of the life-time of (PP)q u - 
decay. There is, however, a serious theoretical problem in the determination of | (m) | from 
experimental data caused by the necessity to calculate the nuclear matrix elements. 

In the framework of Majorana neutrino mixing, the total probability of (PP) o-decay 
has the general form (see Ref. nsa) 

T°"(A,Z) = |{m)| 2 |Af(A,Z)| 2 G ( >‘'(B l) ,Z), (5.26) 

where M(A,Z) is the nuclear matrix element and G 0u (E 0 , Z) is a known phase-space 
factor ( Eq is the energy release). Thus, in order to determine | (m) | from the experimental 
data we need to know the nuclear matrix element M(A, Z). This last quantity must be 
calculated. 

There are at present large uncertainties in the calculations of the nuclear matrix 
elements of (/3/3) 0l ,-decay (see Refs. |152U153U154j ). Two basic approaches to the calcu¬ 
lation are used: quasi-particle random phase approximation and the nuclear shell model. 
Different calculations of the lifetime of the (P/3) oj,-decay differ by about one order of 
magnitude. For example, for the lifetime of the (PP) o^-decay of ,6 Ge, assuming that 
| (m) | = 5 x 10 -2 eV, the range 

6.8 x 10 26 yr < T^^Ge) < 70.8 x 10 26 yr (5.27) 
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A( 76 Ge/ 130 Te) 

11.3 

3 

20 

4.6 

3.6 

4.2 

A( 76 Ge/ 136 Xe) 


1.5 

4.2 

1.1 

0.6 

2 

A( 76 Ge/ 100 Mo) 



14 

1.8 

10.7 

0.9 


Table 2: The results of the calculation of the ratios of the lifetime of (/?/?) 0iy -decay of 
several nuclei in six different models. The references to the corresponding papers are 
given in brackets. 


has been obtained (see Ref. TK31 ). 

The problem of the calculation of the nuclear matrix elements of neutrinoless double-/? 
decay is a real theoretical challenge, ft is obvious that without a solution of this problem 
the effective Majorana neutrino mass | (m) | cannot be determined from the experimental 
data with reliable accuracy (see the discussion in Ref. T55/ T48l| h 

The authors of Ref. [ Il56lj proposed a method which allows to check the results of 
the calculations of the nuclear matrix elements of the (P/3 )oydecay of different nuclei by 
confronting them with experimental data. Let us take into account that 

1. For small neutrino masses (m, < lOMeV) the nuclear matrix elements do not 
depend on the neutrino masses HEp. 

2 . A sensitivity of a few 10~ 2 eV for |(m)| is planned to be reached in future experi¬ 
ments on the search for neutrinoless double-/? decay of different nuclei. 

From Eq. (EH) we have 


R(A, Z/A, Z') 


T°J 2 ( A,Z) 

T°J 2 (A'iZ') 


\M(A',Z , )\ 2 GP v (E , q ,Z') 
\M(A, Z)\ 2 G 0v (E 0 , Z) 


(5.28) 


Thus, if the neutrinoless double /? -decay of different nuclei is observed, the calculated 
ratios of the corresponding squared nuclear matrix elements can be confronted with the 
experimental values. Table El shows the ratios of the lifetimes of (/?/?)oi/-decay of sev¬ 
eral nuclei, calculated in six different models, using the values of the lifetimes given in 
Ref. 11541 . As it is seen from Table El the calculated ratios vary within about one order 
of magnitude. 

As one can see from Table El the values of the ratio i?( 76 Ge/ 130 Te) calculated in 
Ref. |129] and Ref. [Ifilj are, correspondingly, 4.6 and 3.6. It is clear that it will be difficult 
to distinguish models P22! and m through the observation of the neutrinoless double-/? 
decay of 76 Ge and 130 Te. However, it will be possible to distinguish the corresponding 
models through the observation of the (/?/?) 0l/ -decay of ,6 Ge and 100 Mo (the values of 
the corresponding ratio are 1.8 and 10.7, respectively). This example illustrates the 
importance of the investigation of (/?/?)oi/-decay of more than two nuclei. 

The nuclear part of the matrix element of (/?/?)ov-decay is determined by the matrix 
element of the T-product of two hadronic charged currents connected by the propagator 
of a massless boson. This matrix element cannot be connected with the matrix element of 
any observable process. The method proposed in Ref. j!56| is based only on the smallness 
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of neutrino masses and on the factorization of the neutrino and nuclear parts of the matrix 
element of (/3/3)o u - decay. It requires observation of the (f3/3) o^-decay of different nuclei. 

Let us notice that, if the ratio in Eq. (Mb calculated in some model, is in agreement 
with the experimental data, it could only mean that the model is correct up to a possible 
factor, which does not depend on A and Z (and drops out from the ratio (15.281) 1. Such 
factor was found and calculated in Ref. 163!. where in addition to the usual axial and 
vector terms in the nucleon matrix element pseudoscalar and weak magnetic form factors 
were taken into account. It was shown that in the case of light Majorana neutrinos these 
additional terms lead to a universal reduction of the nuclear matrix elements of (/3/3)$ v - 
decay by about 30 %. This reduction, which practically does not depend on the type of 
nucleus, causes a raise of the value of the effective Majorana mass \(m)\ that could be 
obtained from the results of future experiments. 

6 Cosmology 

Perhaps the best example of the fruitful cross-fertilization of high energy physics and 
cosmology is the momentous constraint by Big-Bang Nucleosynthesis (BBN) |32j on the 
number of light neutrino species. Indeed, the number of effective light degrees of free¬ 
dom affects the expansion rate of the Universe; the larger this number, the larger is 
the expansion rate and hence the higher the freeze out temperature of the weak inter¬ 
actions that inter-convert neutrons and protons. Thus, the neutron to proton ratio is 
correspondingly higher and so is the primordial helium yield. These events took place 
when the temperature of the universe was of the order of 1 MeV and therefore it is clear 
that neutrino masses at the 1 eV scale or less play no significant role in primordial light 
element formation. As a consequence, no relevant information on the absolute value of 
light neutrino masses from those early epochs of the history of the universe can be gained. 
This does not mean, however, that cosmology cannot supply interesting information on 
the neutrino mass issue. Fortunately, we can learn about neutrino mass from various cos¬ 
mological and astrophysical instances as different as the Cosmic Microwave Background 
radiation (CMB), the power spectrum in large scale structure (LSS) surveys, and Lyman 
a (Lya) forest studies. We will address these issues in what follows (see also the reviews 
in Refs. [2flll371 :3T?. 3D; THliHlj ). 

6.1 The Gerstein-Zeldovich limit on neutrino masses 

Before entering the issues mentioned explicitly above, let us present the “classical” cos¬ 
mological bound on the sum of the masses of all neutrino species derived by Gerstein 
and Zcldovich (4B.j37.J- Stable light neutrinos (i.e. relativistic at neutrino decoupling) 
are present in the Universe today with an abundance of about 100 neutrinos and an- 
tineutrinos per cm 3 . If they carry mass and this mass is much larger than the present 
CMB temperature (i.e. m u kT CMB ~3x 10 _4 eV, with T CMB ~ 3K), they contribute 
to the known mass density Q m (relative to the critical density p c = 3 Hq/8tiGn, where 
H 0 is the Hubble constant and Gat is the Newton gravitational constant) associated to 
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nonrelativistic matter (mainly dark). The energy density 14 associated to neutrino mass 
can be thus be written as 



( 6 . 1 ) 


where, as usual, the Hubble constant is parameterized as H 0 = lOOh km/s/Mpc. 

Since observationally fl m h 2 ~ 0.15 and < Q m , it follows that Yhi m i < 14 FV. For 
mass degenerated neutrinos this bound implies that rrii < 5 eV for each species. 

6.2 Microwave Background Anisotropies 

The background radiation first detected by Penzias and Wilson in the late fifties follows 
an almost perfect black-body spectrum at the temperature T 0 = 2.728±0.002 K. This ra¬ 
diation is extremely isotropic so that this temperature on the sky is direction independent 
to a precision of 10 -5 , once the Doppler effect due to the peculiar velocity of the Solar 
System is removed. However the Universe is highly inhomogeneous today and this means 
that it should have been sufficiently inhomogeneous in the past in order that structure 
could grow via gravitational instability. Therefore, density inhomogeneities should give 
rise to temperature inhomogeneities in the sky HEE3- For many years such temperature 
fluctuations in the cosmic background radiation have been searched for until they have 
been finally established at the aforementioned minute 1CD 5 level by COBE [122). It is 
customary to expand the temperature fluctuations A T/T 0 in spherical harmonics 



( 6 . 2 ) 


The coefficients cq m are random variables with zero mean (ai m ) = 0 and variance = 

CiSu'dmm' as required by the statistical isotropy of temperature fluctuations. The C\ s 
form the angular power spectrum and this angular power spectrum is conventionally 
shown when presenting CMB results (see Fig. EHt The CMB photons that we now 
record were last scattered at recombination when the universe was about 300, 000 yr old 
and the redshift was z ~ 1100. So, what we get from the temperature fluctuation spec¬ 
trum is essentially a snapshot of the density inhomogeneities at recombination. It is not 
quite a bonafide picture of the anisotropies at that time because in their way to us the 
cosmic photons should have experienced gravitational redshift by changing gravitational 
potentials, the so-called Integrated Sachs-Wolfe effect (ISW) T69j . and eventually, rescat¬ 
tering by ionized gas in interposed clusters of galaxies (Sunyaev-Zeldovich effect [170 :). 
The primary causes of temperature anisotropy, those present at recombination, are three¬ 
fold. There is an intrinsic source associated to the fact that denser spots are hotter and 
hence photons emerging from those denser regions are bluer. But also, photons in denser 
regions will be redshifted as they climb out of their potential wells (Sachs-Wolfe effect 
(SW) nnn). These are competing effects and it depends on the scale under scrutiny 
that one or the other dominates. And finally, the third source of temperature anisotropy 

14 Neutrino masses relate directly to energy density only if the chemical potential of relic neutrinos is 
negligible. It has been shown in j il6'5) that the neutrino chemical potentials of the three species are very 
small. The cosmological limits on neutrino mass that we will discuss in this review comply with this 


fact. 
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Figure 13: A recent plot of the anisotropy spectrum taken from Ref. [ l!T£j . It shows 
BOOMERANG, DASI, MAXIMA, and CBI data (squares and circles). See Ref. [tlRSj 
and references therein for details. 


generation is associated to Doppler shifts arising from the peculiar motion of matter in 
nnderdense regions being attracted towards overdense regions and from which photons 
are last scattered. We collect the three primary sources in the formula: 


AT 



+ <p — n ■ 


v, 


(6.3) 


where </> is the gravitational potential well, v is the peculiar velocity, and n is a unit vector 
pointing in the direction 9, ip. 

To connect the observational CMB anisotropy data with the underlying cosmological 
model and thus have a handle on the different cosmological parameters one has to work 
out the different pieces in the previous equation in terms of the matter density inho¬ 
mogeneities and peculiar velocity of the photon-baryon fluid at recombination. These 
density fluctuations and peculiar velocity, in turn, have to be obtained from the general 
relativistic equations (and/or their newtonian counterparts, when appropriate) to take 
into account their time evolution from given initial (end of inflation) conditions (adia¬ 
batic) up to recombination. Since density perturbations are supposedly small over the 
whole period of interest (up to recombination), one uses linear perturbation theory to 
deal with the problem which then becomes easy to solve. Indeed, a main feature of the 
linearized theory is that by Fourier transforming from r space into k space, the different 
k modes become mutually independent and therefore the corresponding spatial scales 
evolve independently during the linear era of structure formation. Because each k mode 
corresponds to a different spatial size A (A = 2ir/k), a given mode enters the horizon at a 
given epoch. But crossing the Hubble radius is physically relevant since before crossing a 
scale evolves solely under the rule of gravity and only after horizon crossing are also causal 
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effects operative. A sub-horizon sized fluctuation, therefore, experiences both the grav¬ 
itational pull and the pressure gradients of the photon-baryon fluid. Too much gravity 
pull cannot be counteracted by fluid pressure, hence there is a critical size for a perturba¬ 
tion to stand gravity. Beyond that size (so called Jeans size), collapse is unimpeded but 
below the Jeans size collapse can be halted. This Jeans scale is set by the sound speed 
in the primeval plasma because it determines the distance over which a mechanical re¬ 
sponse of pressure forces can propagate over a gravitational free-fall time and thus restore 
hydrodynamical equilibrium in the fluid. For those under-sized perturbations, acoustic 
oscillations set in: compression is followed by rarefaction and back to compression and so 
forth because the in-falling fluid bounces off every time the pressure of the fluid rises to 
the point where it can halt gravitational in-fall and reverse the process from contraction 
to expansion. Since before recombination the pressure of the baryon-photon plasma is 
dominated by photon pressure, the sound velocity is roughly c s ~ c / a / 3 , close to the 
speed of light. So, during the pre-recombination stages, the sound horizon approximately 
matches the Hubble radius and thus scales entering the horizon before recombination 
undergo acoustic oscillations from that moment onwards. Later, when recombination 
takes place and photons are freed from the baryons to which they were previously tightly 
bound, the different modes (scales) are caught at different phases of their oscillation with 
correspondingly different amplitudes of their density perturbations. These compression 
and rarefaction phases translate into peaks in the temperature power spectrum that one 
observes. Odd peaks correspond to compression maxima and even peaks correspond to 
compression minima (rarefaction peaks). The first peak is associated to the scale that 
enters the horizon at recombination and is thus caught in its first oscillation height. The 
second peak corresponds to a scale that has already gone through a complete oscillation 
cycle at recombination, etc. Because the smaller the scale the sooner it entered the hori¬ 
zon, and therefore will have got time for a longer period of oscillations before photon 
decoupling, the corresponding peaks in the power spectrum are progressively attenuated 
as compared to the first compression maximum. The main source of damping, which 
is called Silk damping cm. is due to the fact that photons in the baryon-photon fluid 
have a mean free path governed by the Thomson cross section (photons are coupled to 
the electrons via Thomson scattering and electrons in turn are tightly bound to protons 
via Coulomb interactions) and so photons tend to leak out from overdense regions to 
less dense regions whenever the photon mean free path (which depends on the ionization 
history before recombination and on the baryon content) exceeds the scale of the density 
fluctuation. In addition to this there is also a limit to the pattern of peaks supplied by 
the finite width of the last scattering surface. Since recombination is not instantaneous 
but takes a finite amount of time, observations of the cosmic background temperature 
are actually an average over temperatures of photons that reach us from a shell whose 
thickness is about one tenth the Hubble distance at recombination P22|. Hence, scales 
that are of this order of magnitude or less are completely washed out by the temperature 
averaging process. For a flat universe this limit corresponds to angular scales of about 
0.1 deg or to multipoles larger than about l ~ 2000. We are prepared now to discuss 
what can be learned from the peak structure of the power spectrum as far as general 
cosmological parameters is concerned (including the neutrino energy content, which is 
our main concern here). 

The characteristics of the spectrum of CMB anisotropies and of its peaks ( i.e. their 
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positions, heights, and shapes) depend on the adopted cosmological model for the initial 
anisotropies and on all the parameters that govern the evolution of the acoustic oscilla¬ 
tions before recombination and on the hindrances encountered by photons in their paths 
from recombination onwards. Clearly, not all parameters affect to the same extent the 
different aspects of the power spectrum pattern. Let us remind here of the main influ¬ 
ences on the power spectrum features and let us leave for later the neutrino mass related 
issues: 

1 . At large angular scales (i.e., large compared to the horizon scale at recombina¬ 
tion), the SW effect dominates the anisotropies. Those are the scales that, at 
recombination, have had no microphysics processing yet and reflect directly the 
character and strength of the primordial inhomogeneities. Therefore, one can read 
from the observed spectrum at large angular scales (small /’s) the general cosmo¬ 
logical paradigm about the origin of inhomogeneities. For instance, a flat spectrum 
(as it is indeed the case) in this low l regime suggests an adiabatic scale-invariant 
inflationary-like model of primordial anisotropies. 

2. When a physical feature on the last scattering surface at recombination is projected 
on the sky today, the projection depends on the distance to last scattering and on 
the curvature of space. The distance is fixed by t rec and hence depends mainly on 
the Hubble constant and total energy density Curvature affects the projection 
because it is an obvious geometric fact of curved spaces that objects seem to be 
smaller (larger) in an open (closed) Universe then they would in a flat LIniverse. 
Thus the angle 9 subtended by a given structure on the sky (its angular size) depends 
on curvature. Recall that the first peak corresponds to the scale that enters the 
horizon at recombination. So, its location in l (/ ~ 9~ 1 ) tells us about the curvature 
of the Universe. For a zero curvature Universe the position of the first peak should 
be at l\ ~ 200, which is actually the case of the observed spectrum (see Fig. EEBD- 

3. Baryons contribute inertia to the oscillating photon-baryon plasma in the dark 
matter potential wells and as a consequence the bigger the baryon content is the 
larger the amplitudes of the acoustic vibrations are (i.e., the peaks are higher). 
Furthermore, compression peaks (odd peaks) are enhanced relatively to rarefaction 
peaks (even peaks) as a result of baryon drag. 

Unfortunately there is a substantial amount of degeneracy among the cosmological 
parameters ma that allow for a multiplicity of different parameter choices giving an 
equally acceptable spectrum. So it is very desirable to use alternative measurements 
as complementary tools for determining cosmological parameters and thus help break 
degeneracies. For instance, for a flat LIniverse (f2 ~ 1, as it is indeed the case; see point [21 
above) the position of the first peak is almost independent of the relative weight of matter 
(baryonic plus dark) and dark energy (cosmological constant) in 

Nonzero mass neutrinos affect the CMB anisotropy spectrum to a much lesser extent 
than the previously stated effects scat Their influence is twofold: the position of the 
peaks is slightly modified and also their amplitudes are enhanced. Although the position 
of the first peak is mostly dictated by curvature, the peaks and the troughs move slightly 
to lower I's (to the left in Fig. PHI due to massive neutrinos. This effect can be traced back 
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to the fact that neutrinos being massive, they start being ultrarelativistic until their freeze 
out and beyond and only later in the history of the Universe the neutrinos become nonrel- 
ativistic. Compared to the massless neutrino case where neutrinos are always relativistic 
degrees of freedom, in the case under discussion, the expansion rate is slightly modified 
since the competition between radiation domination and matter domination is altered. 
While the propagation of sound in the photon-baryon fluid depends only on the baryon 
density and hence neutrino mass is not relevant here, the sound horizon at decoupling 
is modified simply because decoupling is slightly delayed due to the change in expansion 
rate. If the sound horizon is larger, so it is the scale that enters the horizon at last scat¬ 
tering. Therefore, larger angular scales corresponding to lower I's, the pattern of acoustic 
peaks is shifted to the left. The neutrino mass, on the other hand, has a larger impact 
on the power spectrum than the shifting just mentioned; it leads to an enhancement of 
the peaks. The origin of the effect is related to the time variation of the gravitational 
potentials. In a radiation dominated Universe potentials change with time whereas in 
a matter dominated Universe gravitational potentials are constant. Since a Universe 
with massive neutrinos implies that relativistic matter turns into non-relativistic matter 
during relevant periods of its evolution, the acoustic oscillations of the baryon-photon 
plasma are being forced by time decaying potentials that differ from those associated to 
a Universe with massless neutrinos only. As a result the acoustic oscillations (mainly 
for l > 300) get an extra boost in amplitude at last scattering (parametric resonance). 
There is also a smaller effect associated to varying potentials after last scattering that 
introduces a relative difference between models with/out massive neutrinos (affecting the 
ISW contribution to anisotropies) which is operative at smaller I's. 

6.3 Galaxy Redshift Surveys 

In the previous section we gave a brief and general description of the physics of the CMB 
angular power spectrum and noted that the direct influence upon it of neutrino mass is 
only marginal. Yet, the CMB power spectrum data is important in the determination 
of neutrino mass because it can be used in combination with other astrophysical sources 
- where the neutrino mass plays a more relevant role - to help reduce the number of 
uncertainties in the various cosmological parameters. One of these sources is large scale 
structure. Neutrino mass affects large scale structure formation and its effect can be 
studied via observation of the distribution of galaxies. Since the distribution of galaxies 
should trace the matter density of the Universe (related to each other via a bias factor), 
large samples of galaxy redshifts in surveys such as the 2 degree Field Galaxy Redshift 
Survey (2dFGRS) (HSJ provide a tool to study the power spectrum of matter fluctuations 
with very small random errors. 

Recall that in linear theory what one is dealing with is the Fourier transformed density 
perturbations 5 k- The initial conditions for 5k are set to reproduce a property of inflation 
(and consistent with observations), i.e. a flat or Harrison-Zeldovich spectrum is assumed. 
This implies that the power spectrum behaves as P(k ) = ( 51) ~ k n with n ~ 1. This 
initial spectrum has to be evolved from the very high initial redshifts to the redshifts 
relevant for structure formation surveys (the median redshift of 2dFGRS is ~ 0.1 C2SD- 
This processing is dictated by the continuity, Euler and Boltzmann equations that govern 
the physics of the perturbations of the cosmic fluid. To be specific, what concerns us 
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Figure 14: 



Ref. [1TE . 


here is the effect of massive neutrinos on the evolution of perturbations, i.e on the power 
spectrum P{k). Once decoupled very early in the history of the Universe (at T ~ 1 MeV) 
neutrinos free-stream at almost the speed of light. This is so until after their momenta 
become on the order of their mass and less and hence they enter a non-relativistic regime. 
During their relativistic life-span they outflow from regions smaller than the horizon so 
that these regions are being depleted and hence energy density perturbations at those 
scales are effectively erased. This phenomenon comes to an end when neutrinos cease to 
free-stream as they become non-relativistic and can cluster with the cold components of 
dark matter for all scales that are larger than the Hubble radius at the time the neutrinos 
become non-relativistic. This limiting scale is given by the formula cm 

k nr ~ 0.03(m i/ /l eV) 1 ^ 2 D^ 2 hMpc _1 . (6.4) 


For all scales smaller than this (i.e. for k > k nr ) the growth of perturbations is suppressed. 
So neutrino mass influences the power spectrum of cosmological structure at small scales. 
The loss of power on small scales can be approximated by mg 


AP <i„ 

- £3 —8 -. 

P o 


(6.5) 


This equation gives us a handle for extracting the bounds on neutrino mass from the 
large samples of data in present and upcoming galaxy distribution surveys (see Fig. ED- 
The analysis of the data has to be restricted to a band of scales for which the data are 
precise enough and for which linear perturbation theory holds. On the large scales side 
accuracy is volume limited and on the small scales side linear theory is jeopardized. The 
2dFGRS data, for instance, are robust on scales 0.02 < k < 0.15 h Mpc -1 . 


6.4 Lyman a forests 

The last piece of information relevant to the neutrino mass issue that we want to discuss is 
the Lycc forest m measurements of the power spectrum of mass fluctuations. Quasars 


43 





















are of help in cosmology because for one thing they rank among the oldest detected 
objects in the Universe and hence provide crucial hints for structure formation studies. 
Furthermore, quasar spectra are a means for studying the intergalactic medium. Atomic 
hydrogen in the gas clouds in the vicinity of the quasar makes 2 p —> Is Lya transitions 
that are seen redshifted today by a factor 1 + z quasar . But what is most important here 
is that quasar spectra show a series of absorption lines associated to intervening clouds 
that photons encounter in the way from the quasar to us. Indeed, the spectrum displays 
a “forest” of Lya absorption lines to the left of the Lya emission line of the quasar, 
i.e. blueshifted relatively to the emitter, that correspond to the resonant absorption 
of those photons - with wavelengths A stretched by cosmic expansion in the proportion 
(1 + z quasar )/(I + z c i ou d ) - that exactly match the Lya transition. The forest of absorption 
lines in the quasar spectrum refers therefore to a sequence of clouds at various redshifts 
along our line-of-sight towards the quasar that absorb radiation from the quasar at specific 
wavelengths given by the redshift of each cloud relative to the quasar’s own redshift. The 
distribution of such clouds thus can provide astronomers with important clues about 
structure formation. In particular, since light neutrinos do not cluster on small scales 
as was explained in the previous section, the measurement of cluster formation on small 
scales extracted from Lya forest observations can lead to definite predictions for neutrino 
mass in the eV range. 

6.5 Neutrino mass bounds 

Now that we discussed the different cosmological and astrophysical sources of information 
on neutrino mass, we can summarize the constraints on neutrino masses that follow from 
these sources ESI- 

Absolute neutrino masses cannot be measured in neutrino oscillation experiments. 
Only mass-squared differences have been established so far. If the three neutrino mass 
spectrum is hierarchical then m i ~ 0, m 2 — \J Arag 0l , and m 3 ~ a/ Am^ tm . Should, 
on the other hand, the mass hierarchy be inverted, then mi ~ 0 , m 2 — a/A rn^m, and 
m 3 ~ ■>/Arn^tm• The third possibility is mass degeneracy, i.e. mi ~ m 2 ~ m 3 . In this case 
the three masses would be much bigger than a /A-rn^ tm and any hint as to their absolute 
mass scale is lost in neutrino oscillations. In this instance astrophysical/cosmological 
tests come to the rescue. 

The Lya forest seen in quasar spectra, as mentioned in the previous section, can be 
used to study mass distribution fluctuations. There is a well-understood theory of Lyct 
forest formation embodied in the standard cosmology which establishes a rather simple 
local connection between the absorbed flux in a quasar spectrum and the underlying 
matter density. From this relationship, the power spectrum P(k ) can be extracted, at 
least over a limited range of scales H3]. 

Indeed, in the usual cosmological scenario where structure formation proceeds via 
gravitational instabilities, the behavior of gas plus a background of UV ionizing radiation 
leads naturally to quasar absorption effects. The physics of the gas is driven by the 
competition of two phenomena, namely adiabatic cooling due to Hubble flow, and heating 
by the photo-ionizing UV background. Hydrodynamic simulations show that the Lya 
forest arises in gas of moderate overdensity 1 tsl) I. This density field can be then locally 
related to the Lya optical depth USD, and consequently to the observed transmitted flux 
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in a quasar spectrum. 

Light neutrinos delay the growth of perturbations at small scales and therefore a 
constraint on their mass is made possible through the recovery of the power spectrum P(k) 
from an observational measurement of the Lya forest. The authors of Ref. used the 
measurement in |18d| to restrict light neutrino masses. Their general strategy consisted of 
two steps. In a first step they tested hydrodynamic Lya forest simulations in the context 
of a cosmological model with non-zero mass neutrinos. Specifically, what they wanted to 
check is their ability to recover, in a given cosmological setting, the power spectrum P(k) 
from a hydrodynamic simulation of the forest. The method had been previously used in 
models that do not include “hot particles”, i.e. models without massive neutrinos. The 
assumed underlying cosmological scenario is an adiabatic, cold dark matter dominated 
Universe with gaussian initial fluctuations. Such models contain six free parameters, 
namely: the matter density fl m , the Hubble parameter h, the baryon density fthe 
neutrino density and the amplitude and tilt of the initial power spectrum of density 
perturbations Ak n . A flat spatial geometry is also assumed throughout as the CMB 
anisotropy and supernova measurements seem to corroborate. The particular model used 
in the tests has Q m — 1, h = 0.5, f \ = 0.075, and 1% = 0.2. The hydrodynamic 
simulation EH3 allows then to follow the evolution of structure in this model and to 
know the physical condition and distribution of the gas at z — 2.5 (i.e. at the redshift 
of the actual measurements). From this, they generated artificial Lya spectra for 1200 
random lines of sight through the simulation volume. The recovery procedure was applied 
then to these spectra to see if the recovered P(k) agrees with the linear theory prediction 
of the model under scrutiny. The authors obtained a fair amount of consistency over the 
whole observational band of scales and concluded from their simulation tests that their 
recovered P(k) was systematically too low in amplitude and had a somewhat flatter slope 
than the linear theory prediction (however, the authors point out that correcting for this 
underestimation would lead only to tighter neutrino mass limits). Having proved that 
the method works for a specific model, it was assumed that it would also work for all 
cosmological models with massive neutrinos available in parameter space. The final step 
involves using the observational results on the Lya forest to explore the six dimensional 
parameter space and place an upper limit on the neutrino mass. 

The loss of power in the Lya forest induced by a neutrino mass is given by equa¬ 
tion (16.51) but other parameters could produce similar power suppression effects. To 
avoid as much as possible undesired degeneracies, the authors use additional cosmologi¬ 
cal constraints. They use the Hubble constant measurement jl85] . the COBE detection 
of large scale anisotropies |186| . the present abundance of galaxy clusters USE, galaxy 
surveys JKS1 , nucleosynthesis limits on baryon abundance [TlROl , and the age of the oldest 
globular clusters to set a lower limit to the age of the Universe USD]. In their analysis 
the authors reject every model that violates the 95% C.L. on any of the aforementioned 
constraints. With such reduced parameter space, the analytic approximations of inn 
were used to find the model that maximizes as a function of Q m . The result is 

< 5.5 eV (95% C.L.), (6.6) 

i 

independently of the value of Q m . 

If Q m is restricted to lie in the range 0.2 — 0.5 as favored by observation, then the 
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authors parameterize their bounds as 15 

m v < 2.4 eV (95% C.L.). (6.7) 

There have been several analysis that put constraints on neutrino masses from LSS 
redshift surveys. The most recent ones make use of the 2dF galaxy survey [175i . The 
2dFGRS is a sample of over 220,000 galaxy redshifts that permits the measurement 
of large-scale structure statistics with very small random errors. Ref. HH2| computes 
the matter power spectrum from linear theory for a multiplicity of cosmological models 
described in terms of the components: baryons, cold dark matter, cosmological constant, 
and of course hot dark matter (i.e. neutrinos with non-zero mass). The calculated matter 
power spectrum and the measured galaxy power spectrum can be put in correspondence 
through a bias parameter, i.e. b 2 = P g [k) / P m (k ), where P g [k ) is the measured galaxy 
distribution power spectrum and P m {k ) is the matter power spectrum. Although b is in 
principle scale dependent, there are good reasons to believe that b is a constant on the 
scales considered in the analysis 193'. The authors of this study absorbed this constant in 
the amplitude A of the power spectrum of density fluctuations taken as a free parameter. 
There is a vast parameter space available for the analysis, and again it helps to take other 
cosmological inputs into consideration. In this way the implications for neutrino mass 
will be less uncertain. 

From primordial nucleosynthesis one has the constraint Vt^ti 2 = 0.020 ± 0.002 on the 
density of baryons (19lj. The Hubble parameter h as measured by the Hubble Space 
Telescope (HST) key project |lffi)] is h = 0.70 ± 0.07. Another prior is the total matter 
density Q m . As stated before when we discussed the peak structure of the CMB anisotropy 
spectrum, there is strong evidence for a spatially flat Universe jl96l . This means Vt m + 
12a = 1. This last relation, used together with the results from surveys of high redshift 
Type la supernovae (SNIa) |197U198j . leads to the constraint = 0.28 ± 0.14. On the 
other hand, independent studies give a wider spread of values. For instance, mass-to- 
light ratio studies of galaxy clusters render typically lower values for Q m (~ 0.15) [TTTTj 
whereas cluster abundance studies deliver Q m ~ 0.3 — 0.9 [ 5UR . Given these facts that 
make Vt m the most poorly known parameter, the authors of the present analysis employed 
two kinds of priors on Vt m . One was a Gaussian at il rn = 0.28 and standard deviation 0.14 
as required by supernova and CMB results and the other was a uniform prior in the range 
0.1 < Q m < 0.5. The latter upper limit (I2 m < 0.5) is dictated by the values of h used 
which would imply, for Q m > 0.5, an age of the Universe shorter than 12 Gyr. Although 
the value n — 1 for the spectral index is the usual theoretical choice, n = 1 ± 0.1 is also 
acceptable and consistent with the CMB data. Therefore, the authors of Reference |l'9'2j 
considered the cases n = 0.9 and n — 1.1 and they ran a grid of models with n as an 
added parameter restricted to the values n — 1 ± 0.1 (Gaussian prior). 

Their results can be summarized as 

J2 m i< 1-8eV (95% C.L.), ( 6 . 8 ) 

i 

15 Here the authors give the parameterization for a single massive neutrino species of mass in the 
degenerate three neutrino case, one should interpret m„ = Jfo rrii. 


46 





for Q m h 2 = 0.15 (for the central values of Vt m and h = 0.7; actually, almost identical 
results follow from the two distinct priors on fl m ) and spectral index n — 1. If this latter 
condition is relaxed to n = 1.0 ± 0.1, then a somewhat looser bound is obtained: 

J2 m i< 2.2 eV (95% C.L.). (6.9) 

i 

Other groups |20112'0'21 obtain results from LSS that are fully compatible with the previous 
upper bounds on neutrino masses. Hannestad in Ref. |M| performs a full numerical 
likelihood analysis in a cosmological parameter space with the following free parameters 
(other than 1%): Q m , h, n, the normalization of the CMB power spectrum Q and 
the optical depth to reionization r. Also, he restricts the study to flat models, i.e. with 
zero curvature. This is by no means a drawback since there is little degeneracy between 
neutrino mass and curvature. The analysis is presented in three stages, depending on 
the data sets used and on the corresponding priors for the cosmological parameters other 
than neutrino mass. In the first stage LSS j2P7>l data and CMB jTP71l2THj data alone are 
used. The resulting bound is, in this instance, 

^mi< 2.96eV (95% C.L.). (6.10) 

i 

The second stage incorporates Big-Bang Nucleosynthesis |205j and Hubble Space Tele¬ 
scope |195| data in the analysis (in addition to the previous data sets). This entails the 
Big-Bang Nucleosynthesis prior on the baryon density fl^h 2 = 0.020 ±0.002 and the HST 
key project prior on the Hubble parameter already given above. This leads to a slight 
improvement of the bound: 

J2 m i< 2.65eV (95% C.L.). (6.11) 

i 

Finally, including the data from high redshift Type la supernova surveys hszi Hannes¬ 
tad obtains: 

^m;<2.47eV (95% C.L.). (6.12) 

i 

Here, the result Q m = 0.28 ± 0.14 (valid for a flat Lhiiverse) has been used. (In all three 
cases above, Q and b are allowed to vary freely, r = 0 — 1 and n = 0.66 — 1.34. In the first 
and second stages, = 0.1 — 1. In stage one, f \h 2 = 0.008 — 0.040 and h = 0.4 — 1.0.) 

Lewis and Bridle m use the sets of observational data on the CMB, LSS, BBN, 
HST and SNIa, that we are already familiar with, to explore the consequences of a non¬ 
zero neutrino mass under somewhat less restricted assumptions than in the analysis just 
discussed (201] . In particular, they consider nine parameter model universes that include 
parameters that account for a “quintessential” equation of state and tensor contributions 
to the power spectrum (allowed in inflationary models). Perhaps the most distinctive fea¬ 
ture of the present analysis is the use of powerful Markov Chain Monte-Carlo techniques 
to perform a fast and efficient exploration of a high dimensional cosmological parameter 
space. As a result of these methods, reference |202| reports 

< 1.5 eV (95% C.L.). (6.13) 
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For the sake of completeness we should include here the result obtained in Ref. [206]. 
namely (for the case of mass degeneracy): 

rrii < 0.9 eV (6.14) 

for h < 0.8, Q m < 0.4, and an age for the Universe in excess of 11.5 Gyr. However, the 
approach of this work is different from what has been discussed here so far. It is based on 
the matching condition for the cosmic mass density fluctuation power at the COBE scale 
and the matter clustering power at the cluster scale. As the authors put it, the advantage 
of using the cluster abundance information is that it refers to the mass function which 
is not affected by any biasing uncertainties (i.e. to what extent galaxies trace the mass 
distribution). 

To summarize all these findings, we can say that the bound 

^2rrii< 3eV (6.15) 

i 

should be a reliable upper limit on the sum of neutrino masses which implies that (again, 
for three almost degenerate neutrinos) 

nii < 1 eV (6.16) 


for each of the three masses. 

As to future prospects, MAP/PLANCK CMB data in conjunction with high precision 
galaxy surveys such as the Sloan Digital Sky Survey mu could render C23 


or 


y: rrii < 0.3 eV 


< 0.12 eV 


(6.17) 


(6.18) 


at 95% C.L., or even go down to an ultimate sensitivity of about 0.04 eV when weak 
lensing of galaxies by large scale structure is also taken into account 2TT7J . 

To conclude, perhaps we should mention here a result |210| that has an extra the¬ 
oretical input, namely leptogenesis as the origin of matter-antimatter asymmetry. In 
such scenario, neutrinos are Majorana particles, and for the whole picture to work, such 
neutrinos should weigh less than 0.2 eV. 


7 Cosmic Rays 

Since the sixties mm it is well-known that the universe is opaque to protons (and other 
nuclei) on cosmological distances. An ultra high energy (UHE) proton with energy E 
exceeding the Greisen-Zatsepin-Kuzmin (GZK) energy 

E G zk ~5x 10 19 eV (7.1) 

interacts with the photons of the cosmic background producing pions through the A* 
resonance, p + 7 cb —► A* —> Nn. I 11 this way, the initial proton energy is degraded with 
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Figure 15: UHE Cosmic Ray Flux times E 3 . Results from the HiRes-I and HiRes-II 
detectors, and the AGASA experiment are shown. Also shown (solid curve) is a fit to 
the data assuming a model with two sources of cosmic rays, galactic and extragalactic, 
which includes the GZK cutoff. Figure taken from Ref. | 222 j. 

an attenuation length of about 50 Mpc HUSH- The UHE photons have even shorter 
absorption lengths (~ 10 Mpc for E ~ 10 2 ° eV mu) due to their interactions with cosmic 
background photons 2T2 . Since plausible astrophysical sources for UHE particles (like 
AGNs) are located at distances larger than 50-100 Mpc, one expects the so-called GZK 
cutoff in the cosmic ray flux at the energy given by dm 

However, in the cosmic ray data ITHl'2141121511210112171121811219112201221112221 . there are 
about twenty cosmic ray events with energies just above the GZK energy (see Fig. d. 
Yet, the whole observational status in the UHE regime is controversial. While the HiRes 
collaboration claim [ 222 ] that they see the expected event reduction, a recent reevaluation 
of AGASA data seems to confirm the violation of the GZK cutoff |22I]. Indeed, some 
apparent inconsistencies among data have been pointed out 52)1 • The observational 
status is not settled, but it is clear that if the GZK violation is confirmed, the origin of 
the super-GZK particles constitutes one of the most pressing puzzles in modern high- 
energy astrophysics (for a recent review, see Ref. ? 224j ). 

Several hypothetical explanations have been put forward to account for this phe¬ 
nomenon. For example, there are scenarios where the UHE cosmic rays are decay prod¬ 
ucts of exotic super-massive particles or relic topological defects. Also, a way to solve 
the problem might be to postulate a violation of Lorentz symmetry, or introduce new 
particles or/and new interactions (for a recent review, see Ref. [ 224j ). We are concerned 
here with a possible explanation that is based on the so-called Y-bursts [HUE 2 ]- If the 
GZK cutoff is violated and the Z-burst mechanism is indeed the solution to the GZK 
puzzle, it may be used to determine the absolute value of neutrino masses and, in fact, 
it would be an indirect proof of the existence of the relic cosmic background neutrinos. 

The main hypothesis of the Z-burst explanation of the GZK puzzle is the existence of 
a very high flux of UHE neutrinos. And the main criticism to it is that standard astro- 
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physical objects cannot produce such fluxes. Thus, the Z-burst hypothesis requires new 
sources producing UHE neutrinos copiously. (If the flux is so large, it may be measured 
by the next generation of neutrino telescopes; see Ref. |225| .) However, compared to the 
other models quoted before, the Z-burst scenario does not need new physics beyond the 
Standard Model of particle physics with neutrino masses. 

The crucial observation is that neutrinos, contrary to protons and photons, propagate 
over cosmological distances with negligible opacity. The most important interaction they 
have is with cosmic background neutrinos. Although this represents a extremely small 
opacity it might be enough to generate the UHE cosmic rays. 

Cosmic background neutrinos are the relics of neutrino decoupling in the early Uni¬ 
verse, that happened when the Universe was about 1 second old and had a temperature 
of about 1 MeV. Their number density today can be calculated easily in terms of the ob¬ 
served number density of the relic microwave background n 1 (see, for example, Ref. |22B])- 
In the case of no net leptonic number it is given by 


n v = n D . = — n~ = 56 cm 3 . 

u i u i ' ' 


(7.2) 


where z/j is a neutrino species. The neutrino cosmic background is in many aspects similar 
to the photon cosmic background. The photons have been detected but the neutrinos not 
yet, because the photons have much stronger interactions than the neutrinos. In any case, 
the neutrino background is a firm prediction of cosmology. These neutrinos decoupled 
when relativistic, but today they are non-relativistic (if the neutrino mass rrq > 1CU 4 
eV.) 

The interaction of UHE neutrinos with background neutrinos is strongly enhanced 
when it proceeds through the Z-resonance, vv —> Z —► hadrons. To be exactly on 
top of the resonance, the UHE neutrino has to have an energy E R such that Mf = 
(E r + mi) 2 — ( E r ) 2 ~ 2miE R : so that 



Ml 

2 mi 


~ 4.2 



x 10 21 eV. 


(7.3) 


We see there is a resonant energy for each neutrino mass 777 , and for this reason we have 
added the subscript z/; to E R . The lowest value of E R corresponds to the highest neutrino 
mass. 

To have the Z-resonance enhancement, the neutrino energy E must be in the range 


E ~ 


£. r ± r. 


(7.4) 


where F z is the Z width. 

The idea of the Z-burst mechanism CUES is that when a UHE neutrino with energy 
E in the band (Q scatters off a relic neutrino, about 70% of the times gives 

vv —> Z —► hadrons . (7.5) 

The hadrons form a highly collimated final state, since in the relic neutrino rest frame 
the Z particle has a Lorentz factor of 7 ~ 10 11 . Provided the scattering takes place at a 
distance from Earth of less than 50 Mpc, attenuation is small and the produced particles 


50 



may induce air showers in the Earth atmosphere giving rise to the observed super-GZK 
events. 

The properties of the process e _ e + —> Z —» hadrons, have been studied with huge 
statistics at LEP (see Ref. 22/ for a review), and can be used to determine the properties 
of the reaction m- The average multiplicity is ~ 30, with final states having an average 
of about 2 nucleons, 10 neutral pions and 17 charged pions. The nucleons will then have 
an average energy 

(E p ) ~ ~ 1.4 ) x 10 20 eV. (7.6) 

30 \rriij 

Neutral pions originate photons, tt° —> 77 , with average energy 

(E~) ~ ^ ~ 0.7 ( x 10 20 eV . (7.7) 

x 7/ 60 \mj v ; 

In the Z-burst scenario these protons and/or photons are the cosmic ray primaries. We 
see that one then needs a mass 777 ~ 0.1 — 1 eV to be just above the GZK cutoff (Hi- 

Relevant for the purposes of this review is the determination of the neutrino mass 
if the Z-burst mechanism turns out to be the solution to the GZK puzzle, as discussed 
in Ref. [2BJ (see also Ref. [2281 ). A detailed analysis has been done in Ref. [ 2223 , where 
the observed UHE cosmic ray spectrum is compared with the predictions of the Z-burst 
model. In the analysis one uses collider data to derive the spectra of the final state, and 
finally one determines the energy losses in the propagation of particles until reaching the 
Earth atmosphere. A maximum likelihood analysis gives the interval 

0.08 eV < mi < 1.3 eV (68% C.L.), (7.8) 

assuming that the super-GZK events are originated by protons produced in Z-bursts 
outside our galaxy. The authors of Ref. [ 2211 ] claim that this neutrino mass determination 
is fairly robust against variations of presently unknown quantities. 

Another interesting study concerns the Z-burst model in presence of a leptonic asym¬ 
metry [ 230j . where one has n Vi % n Pi and (17.21) is no longer valid. The authors of Ref. |230| 
conclude that a neutrino mass m* ~ 0.07 eV, consistent with Super-Kamiokande data, 
explains the cosmic ray events in this leptonic asymmetric case. Also, we would like to 
mention that the possibility that Z-bursts may account for events just below the GZK 
energy and above the ankle of the cosmic ray spectrum is considered in Ref. BHU , again 
for a mass m, ~ 0.07 eV. 

Many other aspects of the Z-burst scenario have been treated in the literature [232]. 
For example, to what extent may neutrino clustering (which is quite likely) enhance 
the signal, which other observations may put constrains to the model, and which are 
the distinctive features that may help us in discriminating between Z-bursts and other 
explanations of the UHE cosmic ray puzzle. 

The future experimental projects aimed at the detection of cosmic rays to probe the 
UHE regime will be crucial to shed more light on this subject (for a recent review, see 
for example Ref. I2H3 T. 


51 


8 Supernova Neutrinos 

Supernovae are extremely powerful explosions that terminate the life of some stars. Typ¬ 
ically some solar masses are ejected in the interstellar space with a kinetic energy of the 
order of 10 51 ergs. The turbulence produced in the stellar medium can help the formation 
of new stars. The ejecta contain heavy elements that are important for the chemical evo¬ 
lution of galaxies, stars, planets and life. Some supernovae produce a compact remnant, a 
neutron star or a black hole, that may be observed. Reviews and extensive bibliographies 
on the physics of supernovae can be found in Refs. [E31I2M4U17] . A fairly updated general 
introduction and review of supernova neutrino physics is given in Refs. Hung. 

Several supernovae that have exploded in parts of our galaxy not obscured by dust 
have been observed with naked eye during the last 2000 years. Most famous is the 1054 
supernova that produced the Crab nebula and the Crab pulsar. The 1006 supernova is 
the brightest supernova of all times. The last galactic supernovae have been observed 
by naked eye in 1572 (Tycho Brahe) and 1604 (Joannes Kepler). In the last centuries 
many supernovae occurring in other galaxies have been observed with telescopes because 
their luminosity is comparable to that of an entire galaxy. Supernova SN1987A, which 
occurred on 23 February 1987 in the Large Magellanic Cloud, is the best studied of all 
supernovae and it is the only one which has been detected also through its neutrino burst. 
As we will see in the following, this first historical observation of neutrinos produced out 
of the solar system (and even out of our galaxy) is important not only for the study 
of supernova dynamics, but also for the study of neutrino properties, and in particular 
neutrino mass. 

8.1 Supernova Types and Rates 

For historical reasons, supernovae are divided in the four different types listed in Table 01 
characterized by their spectroscopic characteristics near maximum luminosity, which de¬ 
pend on the composition of the envelope of the supernova progenitor star. The two wide 
categories called Type I and Type II are characterized by the absence or presence of hydro¬ 
gen. However, the most important physical characteristic is the mechanism that generate 
the supernova, that distinguishes supernovae of Type la from supernovae of Type lb, Ic 
and II, as shown in Table 01 This difference becomes noticeable from the electromagnetic 
spectrum some months after maximum luminosity, when the ejecta become optically thin 
and the innermost regions become visible. 

Typically the optical emission of both Type I and II supernovae start with a rise in 
luminosity during a week or two, due to the expansion of the luminous surface. Type I 
supernovae have typically a narrow luminosity peak, whereas Type II have broad peaks, 
of the order of 100 days. After the peak the luminosity decreases during about one year. 

Type la supernovae are thought to be generated by carbon-oxygen white dwarfs that 
have a close companion star from which the white dwarf can accrete mass. When the 
mass of the white dwarf reaches the Chandrasekhar limit 16 of about 1.4M 0 , the star 

16 White dwarfs are the evolutionary product of stars that have finished thermonuclear fuel burning. 
They weight about one solar mass, they have a radius of about 5000 km and a density of the order of 
10 6 g cm -3 . The pressure of degenerate electrons support white dwarfs against the inward pull of gravity 
(see, for example, Ref. 12351 b In 1931 Chandrasekhar discovered that white dwarfs have a maximum 
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Type 

near maximum 

months later 

Mechanism 

Remnant 

H 

He 

Si 

Fe 

0 and C 

la 

No 

No 

Yes 

Yes 

No 

Mass 

Accretion 

None 

lb 

Yes 

No 

No 

Yes 

Core 

Collapse 

Neutron Star 

or 

Black Hole 

Ic 

No 

No 

II 

Yes 

? 

? 


Table 3: Main characteristics of supernova types. 

becomes unstable, because the pressure of the degenerate electron gas cannot sustain any 
more the gravitational weight. The white dwarf begins to collapse, triggering the fusion of 
carbon and oxygen to heavy nuclei, that liberate an enormous quantity of energy causing 
the explosion of the star (see Ref. |23dj h This explosion disrupts the progenitor white 
dwarf and generates an expanding nebula without a central compact object. 

Since supernovae of Type la are all generated under similar physical circumstances, 
they have almost identical characteristics, the most important being the total luminosity 
and the “light curve” (luminosity as a function of time). An empirical relation between the 
duration of the peak phase of the light curve and the luminosity of Type la supernovae has 
been discovered by Phillips in 1993 [2H7J from the catalog of observed Type la supernovae 
in nearby galaxies with known distance. This width-luminosity relation (“broader is 
brighter”) allows to use Type la supernovae as standard candles for the measurement of 
the distance of galaxies as far as 100 Mpc or more (see Ref. |2HBj and references therein). 

The observation by the Hubble Space Telescope of supernovae of Type la in galaxies 
at cosmological distances have recently been used for the measurement of the Hubble 
parameter and the deceleration constant. Contrary to the expectations, it has been 
found that the rate of expansion of the Universe is accelerating {23911240] . This surprising 
behavior can be explained in the framework of the standard Friedmann-Robertson-Walker 
cosmology (see Ref. [ 241] ) through the presence of a relatively large vacuum energy (“dark 
energy” or a cosmological constant). 

From the point of view of neutrino physics, Type lb, Ic and II supernovae are much 
more interesting than Type la supernovae, because they produce a huge flux of neutrinos 
of all types. This is due to the fact that these supernovae originate from the collapse of 
the core of massive stars (M > 8M 0 ) that leaves a compact remnant. During the few 
seconds following the collapse, the compact remnant is very hot and neutrinos of all types 
are copiously produced. Since the remnant and the surrounding envelope are optically 
thick, about 99% of the gravitational binding energy liberated by the collapse (about 
3 x 10 53 ergs) is carried away by neutrinos. The average energy of the emitted neutrinos 
and antineutrinos is of the order of 10 MeV, and their number is about 10 58 , about one 

mass of about 1.4 Mq, above which the star collapses. 
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galaxy 

supernova type 

type 

la 

lb, Ic 

II 

All 

E- SO 

0.32 ±0.11 

< 0.02 

< 0.04 

0.32 ±0.11 

SOa - Sb 

0.32 ±0.12 

0.20 ±0.11 

0.75 ±0.34 

1.28 ±0.37 

Sbc - Sd 

0.37 ±0.14 

0.25 ±0.12 

1.53 ±0.62 

2.15 ±0.66 

All 

0.36 ±0.11 

0.14 ±0.07 

0.71 ±0.34 

1.21 ±0.36 


Table 4: Supernova rates in units of h 2 SNu from Refs. [24411242] . One supernova unit 
(SNu) is defined as one supernova per 10 10 per 100 yr, where T 0 ,b is the solar 
luminosity in the blue spectral band. 


order of magnitude larger than the lepton number of the collapsed core. 

Type 11 supernovae are though to be generated by the core collapse of red (or blue 
as SN1987A) giant stars with a mass between about 8 and 60 solar masses. Since the 
size and mass of the hydrogen envelope can be very different from star to star, even 
if they have the same initial mass, the visible effects of the supernova explosion have 
a wide range of variability, leading to a further classification of Type If supernovae as 
Type 11L if the decrease of the luminosity is approximately linear in time, as Type IIP 
if the time evolution of the luminosity shows a plateau, as Type IIF if the supernova is 
faint, as Type lib if helium dominates over hydrogen, as Type Iln if the spectrum shows 
narrow line emissions, as Type Ilpec if the supernova has peculiar characteristics (see 
Refs. |242ll243j : subclasses determined by spectral properties are denoted by lower-case 
letters and subclasses determined by properties of the light curve are denoted by upper¬ 
case letters). It is believed that if the exploding star does not have a hydrogen envelope 
the supernova is of Type lb, and if also the helium shell is missing the supernova is of 
Type Ic. All these classes are not clear-cut and intermediate cases exist. 

Supernova SN1987A was an extreme case of Type IIP, since the luminosity increased 
for about 3 months after collapse and the supernova was rather faint. Therefore, some¬ 
times SN1987A is classified as IIP |242U243| . sometimes as IIF [ 235 and sometimes as 
Ilpec |243j . It is believed that its faintness is due to the compactness of the progenitor 
(a radius of about 10 12 cm). In this case much of the available energy is used in the ex¬ 
pansion and the luminosity increases for some time because of the growing contribution 
of radioactive decay of heavy elements in inner shells, that become more visible as the 
envelope expands. 

A very important problem is the estimation of supernova rates. Table Q] shows the 
recent estimates of supernova rates presented in Refs. 123212323, that have been obtained 
from the Asiago Supernova Catalog 123512351 - Some of these rates are significantly smaller 
than previously thought [2371 . One can see that the rate of core-collapse supernovae of 
Type lb, Ic and II depends rather strongly on the galaxy type, being very small in elliptical 
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phase 

1 Mq 

25 M q 

(£) 

p c 

(g cm -3 ) 

At 

(yr) 


Pc 

(g cm~ 3 ) 

At 

(y r ) 

H burning 

1.3 

153 

1.1 x 10 1U 

3.3 

3.8 

6.7 x 10 6 

He burning 

11 

2.0 x 10 4 

1.1 x 10 8 

17 

762 

8.4 x 10 5 

C burning 




72 

1.3 x 10 5 

522 

Ne burning 




135 

4.0 x 10 6 

0.89 

O burning 




180 

3.6 x 10 6 

0.40 

Si burning 




314 

3.0 x 10 7 

2.0 x 10“ 3 


Tabic 5: Central temperature T c , central density p c and time scale At of the evolutionary 
phases of Population 1 stars with initial masses 1M 0 and 25 M 0 (values taken from 
Ref. 20) )• 


galaxies. These galaxies are very old and have little star formation that could produce 
short-lived massive stars that end their life with a core-collapse supernova explosion. 
Instead, the rate of Type la supernovae is almost independent from the galaxy type, 
because mass accretion can occur also in old population II star. 

One of the most crucial questions for supernova neutrino astronomy is the rate of 
core-collapse supernovae in our galaxy, that could produce an observable neutrino burst 
with high statistics in neutrino telescopes. The morphological type of the Milky Way 
is thought to be Sb-Sbc and the luminosity is 2.3 x 10 10 L Q) b- From Table U using 
a Hubble parameter h ~ 0.7, the rate of core-collapse supernovae in the Milky Way is 
about 2±1 per century. This rate is about a factor of two smaller than previous estimates 
derived from counts of historical supernovae and of supernovae remnants |247j . but the 
large uncertainties do not allow to claim a disagreement and leave the problem open to 
further study. The lack of observation of neutrinos from core-collapse supernova in our 
galaxy since the Baksan Underground Scintillator Telescope began observations in June 
1980 is consistent with the estimated rate and implies that the true rate cannot be much 
higher 12481 . 

8.2 Core-Collapse Supernova Dynamics 

Since only supernovae produced by the collapse of the core of massive stars produce large 
fluxes of neutrinos that could be detected on Earth, here we present a short description of 
the current standard theory of the dynamics of core-collapse supernovae and the resulting 
neutrino flux (see Refs. 1B351250125112521541234117] and references therein). As explained 
in the previous subsection, core-collapse supernovae are classified as of Types II, lb or Ic 
depending on their spectroscopic characteristics at maximum luminosity. However, these 
characteristics depend only on the composition of the envelope, which play no role in the 
collapse of the core and neutrino production. Hence, the following theory applies equally 
well to all Types II, lb, Ic core-collapse supernovae. 

It is believed that core-collapse supernovae are the final stage of the evolution of stars 
with mass between about 8 and 60 solar masses. Lighter stars end their life as white dwarfs 
(but may explode as Type la supernovae if they belong to a multiple system), whereas 
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heavier stars are unstable and probably collapse into black holes without a supernova 
explosion. Stars with mass in excess of 12 solar masses are thought to undergo all the 
stages of nuclear fusion of hydrogen, helium, carbon, neon, oxygen, silicon (see Tableland 
Ref. i 234j h until the star has an onion-like structure shown in Fig. 1 1 01 with an iron core 
surrounded by shells composed by elements with decreasing atomic mass. At this point 
the iron core has a mass of about 1 solar mass, a radius of a few thousand km, a central 
density of about 10 lo gcm -3 , a central temperature of about 1 MeV, and its weight is 
sustained by the pressure of degenerate relativistic electrons. Since iron is the most bound 
nucleus, there is not any more thermonuclear fuel to burn: the iron core is endothermic; 
it can only absorb energy by breaking into lighter nuclei or creating heavier elements. 
The equilibrium between the inward pull of gravity and the electron pressure that sustain 
the core is destabilized shortly before the core has reached the standard Chandrasekhar 
mass of about 1.4 M 0 , because the core contracts and the increased temperature causes 
photodissociation of iron through the process 

7 + 56 Fe —► 13 a + An . (8.1) 

This reaction absorbs about 124 MeV of energy and reduces the kinetic energy and 
pressure of the electrons. Electron capture of nuclei, 

e“ + Af(Z, A) —> Af(Z — 1, A) + z/ e , (8.2) 


and free protons, 

e~ + p —> n + u e , (8.3) 

favored by the high electron Fermi energy, additionally reduces the number and pressure 
of the electrons. At the onset of collapse, when the density of the iron core is not 
too high, the electron neutrinos produced by electron capture leave the core carrying 
away most of the kinetic energy of the captured electrons. The combined effect of iron 
photodissociation and electron capture lowers the electron pressure, decreasing the value 
of the Chandrasekhar mass, until the Chandrasekhar mass becomes smaller than the core 
mass. At this moment the pressure of degenerate relativistic electrons cannot sustain the 
weight of the core any more and collapse commences. As the density and temperature 
increase the processes (EU-flES) proceed faster, lowering further the electron pressure 
and favoring the collapse, which accelerates. 

According to theory (see Ref. 0 and references therein), stars with mass between 
about 8 and 12 solar masses burn hydrogen, helium, carbon, but the core does not get 
hot enough to burn oxygen. However, the core contains neon and magnesium at high 
density, which can undergo electron capture, reducing the electron pressure that sustains 
the core against gravity. As a result, the core collapses and during the collapse oxygen, 
neon and magnesium are converted to iron. Therefore, also in this case the supernova 
explosion is produced by the gravitational energy released by the collapse of an iron core. 

The collapse of the core produces a neutron star and the huge liberated gravitational 
energy is released mainly as a flux of neutrinos, with a small fraction as electromagnetic 
radiation and kinetic energy of the ejecta, which constitute the visible explosion. The 
liberated gravitational energy is about 
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Figure 16: Onion-like interior structure of a Population I star of 25 M 0 just before the 
onset of collapse (see Ref. [ 21151 ). Fe represents assorted iron-peak elements: 48 Ca, 50 Ti, 
54 Fe, 56 Fe, 58 Fe, 66 Ni. The Si shell contains less abundant amounts of S, O, Ar, Ca, 
the O shell contains less abundant amounts of Ne, C, Mg, Si, the He shell contains less 
abundant amounts of C, Ne, O, and the H shell contains less abundant amounts of He, 
Ne, O, N, C. 

where G is Newton gravitational constant, M core is the mass of the core and R is the 
radius of the neutron star. Since M core ~ 1M 0 and R ~ 10 km, the released gravitational 
energy is of the order of 3 x 10 53 ergs. Only about 0.01% of this energy is transformed 
into electromagnetic radiation and about 1% is transformed into kinetic energy of the 
ejecta. 

Let us examine in more detail the mechanism of formation of the neutron star, of 
neutrino production and of supernova explosion. 

The electron neutrinos produced by the electron capture processes m and Q) 
initially leave freely the core, carrying away energy and lepton number, since their mean 
free path is longer than the radius of the core. In this so-called “capture” phase, electron 
neutrinos have a non-thermal spectrum and average energy that grows from about 12 to 
about 16 MeV (see Ref. ; 2K0j ). The luminosity reaches about 10 53 ergssec -1 , but in total 
only about 10 51 ergs are released before core bounce, because the capture phase is too 
short (less than about 10 ms). 

When the density of the inner part of the core (about 0.8 M 0 ) exceeds about 3 x 
10 11 g cm -3 neutrinos are trapped in the collapsing material leading to an adiabatic col¬ 
lapse with constant lepton number. During this stage, the inner part of the core collapses 
homologously, i.e. with subsonic velocity proportional to radius. The outer part of the 
core collapses with supersonic free-fall velocity. 

After about one second from the start of instability, the density of the inner core 
reaches the density of nuclear matter, about 10 14 gcm -3 , and the pressure of degenerate 
non-relativistic nucleons abruptly stops the collapse. The inner core settles into hydro¬ 
static equilibrium, forming a proto-neutron star with a radius of about 10 km, while a 
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supersonic shock wave caused by the halting and rebound of the inner core forms at its 
surface. The shock propagates outward through the outer iron core, which is still col¬ 
lapsing, with an initial velocity of the order of 100 km msec -1 . The gas that is infalling 
at a velocity near free-fall is abruptly decelerated within the shock. Below the shock it 
falls much more slowly on the surface of the proto-neutron star, accreting it. Therefore, 
the proto-neutron star develops an unshocked core with nuclear density, of the order of 
10 14 g cm -3 and radius of the order of 10 km, and a shocked mantle with decreasing den¬ 
sity, down to about 10 9 gcm -3 and a radius of about 100 km, up to the surface of the 
proto-neutron star, where the density has a steep decrease of several orders of magnitude. 

As the shock propagates through the infalling dense matter of the outer core, its energy 
is dissipated by the photodissociation of nuclei into protons and neutrons. Thus, the 
material behind the shock wave is mainly composed of free nucleons. Free protons have a 
high electron capture rate, leading to the transformation of most protons in neutrons, with 
huge production of electron neutrinos. These neutrinos pile up behind the shock, which is 
dense and opaque to them, until the shock reaches a zone with density about 10 11 g cm -3 
(“shock breakout”) a few milliseconds after bounce and the electron neutrinos behind 
the shock are released in a few milliseconds. This neutrino emission is usually called 
“prompt electron neutrino burst” or “neutronization burst”, to be distinguished from the 
thermal production of all neutrino flavors. The neutronization burst has a luminosity of 
about 6 x 10 53 ergs sec -1 and carries away an energy of the order of 10 51 ergs in a few 
milliseconds. In spite of his name, the neutronization burst is too short to carry away a 
significant part of the electron lepton number of the core, which remains trapped. Only 
the low-density periphery of the proto-neutron star is neutronized. 

The energy lost by photodissociation of nuclei and neutrino emission weakens the 
shock (about 1.5 x 10 51 ergs are dissipated for each 0.1 solar masses of photodissociated 
material). In the so-called “prompt” supernova explosion scenario, the shock, although 
somewhat weakened, is able to expel the envelope of the star generating the supernova 
explosion on a time scale of the order of 100 msec. If the star weights more than about 10 
solar masses, the shock is weakened and stalls about 100 ms after bounce, at a radius of 
about 200-300 km, with insufficient energy to reach the outer layers of the star. Matter 
continues to fall through the stalled shock and be photodissociated. If too much matter 
lands on the proto-neutron star, the pressure of degenerate nucleons is not sufficient 
to maintain stability and the core collapses into a black hole, presumably without a 
supernova explosion. 

The conditions that lead to a prompt supernova explosion, without a stalling shock, 
are controversial and are thought to depend on the mass of the progenitor star and on the 
equation of state of nuclear matter, which determines the energy transferred to the shock 
wave by the bounce. It is widely believed that in order to obtain a supernova explosion 
if the shock stalls, the shock must be revived by some mechanism that is able to renew 
its energy. The mechanism which is currently thought to be able to revive the shock is 
the energy deposition by the huge neutrino flux produced thermally in the proto-neutron 
star mm- I n this case, a so-called “delayed” supernova explosion is produced on a 
time scale of the order of 0.5 sec after bounce. 

Neutrinos of all flavors are produced in the hot core of the proto-neutron star (see 
Refs. ). which has a temperature of about 40 MeV, through electron-positron pair 


annihilation, 


e + e + —> v + v , 

(8.5) 

electron-nucleon bremsstrahlung, 


e ± + N^e ± + N + u + P, 

( 8 . 6 ) 

nucleon-nucleon bremsstrahlung, 


N + N —> N + N + v + is , 

(8.7) 

plasmon decay 

7 —> v + 9, 

( 8 . 8 ) 

and photoannihilation 

y + e 1 * 1 —>e ± + z/ + P. 

(8.9) 


Electron neutrinos are also produced by the electron capture process (E3D, and electron 
antineutrinos are produced by positron capture on neutrons (e + + n —> p + 9 e ). In spite 
of their weak interactions, these neutrinos are trapped in a supernova core because of the 
very high matter density. Neutrino can free-stream out of the mantle proto-neutron star 
only at a distance from the center where the matter density is low enough (of the order 
of 10 n gcm~ 3 ) that the mean neutrino free path is larger than the radius of the core. 
The sphere from which neutrinos free stream is called neutrino sphere, and it lies within 
the mantle of the proto-neutron star. Since neutrino interactions depend on flavor and 
energy (see Ref. |255| ). there are different energy-dependent neutrinospheres for different 
flavor neutrinos. More precisely, since the medium is composed by protons, neutrons 
and electrons, and the neutrino energy does not allow creation of muons and taus, the 
flavor neutrinos v e and v e can interact with the medium through both charged-current 
and neutral-current weak processes, whereas the neutrinos u T , v r can interact only 

through neutral-current weak processes, which are flavor-independent (small differences 
between neutrino and antineutrino interactions can be neglected). Therefore, there are 
three energy-dependent neutrinospheres: one for z/ e , one for v e and one for u^, u T) 
v T . From now on we will denote u T , P T collectively as u x , as usually done in 

the literature. Each energy-dependent neutrinosphere emits a black-body thermal flux 
of neutrinos at the considered energy. The estimated radii of the neutrinospheres lie 
between about 50 and 100 km. As we have seen above, when the shock passes through 
the electron neutrino neutrinosphere (shock breakout) a few milliseconds after bounce, a 
large flux of electron neutrinos is released in a few milliseconds in the neutronization burst. 
After shock breakout each neutrinosphere produces a thermal flux of the corresponding 
neutrino flavor. 

The opacities of v e and P e are dominated, respectively, by the charged-current weak 
interaction processes 


v e + n —> p + e , ( 8 . 10 ) 

P e + p —> n + e + . (8.11) 

These reactions allow exchange of energy and lepton number between the neutrinos and 
the medium (which is composed by electrons, positrons, nucleons and photons). For 
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example, in the process (IHTol) the neutrino energy is mainly transferred to the final 
electron 1 ' whose creation increases by one unit the lepton number of the medium. 

Since the mantle of the proto-neutron star is neutron-rich, the opacity of v e of a 
given energy is larger than the opacity of v e with the same energy, and the corresponding 
v e neutrinosphere has larger radius than the v e neutrinosphere. Therefore, for a fixed 
neutrino energy z/ e ’s are emitted by a deeper and hotter layer than z/ e ’s, leading to a 
D e mean energy larger than the v e mean energy. Moreover, the spectra do not have a 
perfect black-body shape, but are “pinched”, i.e. both the low- and high-energy tail are 
suppressed with respect to the tails of a black-body thermal spectrum with the same mean 
energy. Figure El shows the time evolution of neutrino luminosity and average energy 
obtained with the numerical supernova model in Ref. |256j . Other similar estimations 
of the neutrino luminosity and average energy have been obtained with the numerical 
simulations in Refs. [25711258] . A rough estimate of the time-integrated average energies 
is 

(E Ue ) fa lOMeV, (E Pe ) fa 15MeV, (E Ux ) fa 20MeV. (8.12) 

In the delayed supernova explosion scenario the stalled shock lies at a radius of about 
100-300 km, well outside of the neutrinosphere. The post-shock temperature is about 
1.5 MeV and the density of the order of 10 8 gcm -3 . The capture of a small fraction, 
about 5-10% [250] . of the thermal flux of neutrinos emitted from the neutrinosphere 
could revive the shock, leading to the explosion. The largest energy deposition is due to 
electron neutrinos and antineutrinos, which have a charged-current cross section on the 
free nucleons behind the shock that is larger than the neutral-current cross section of all 
neutrino types and is able to deposit more energy. 

If enough energy is deposited behind the shock, about half second after bounce the 
shock is revived and starts to sweep the outer layers of the star generating the explosion. 
Unfortunately, most one-dimensional (i.e. spherically symmetric) computer simulations 
[260112571 1267 did not obtain a successful explosion, which was recently obtained only 
by the Livermore group [ 256) (they used the so-called “neutron finger convection” in 
the proto-neutron star to enhance the early neutrino luminosity which leads to a large 
energy deposition behind the shock). In recent years several groups have performed two- 
dimensional simulations (i.e. cylindrically symmetric) with unsatisfactory results (see 
Ref. [262.. 259j ) and recently a successful three-dimensional simulation of explosion has 
been presented in Ref. |2E3|. The multi-dimensionality of the simulations is important in 
order to take into account convection effects that enhance the efficiency of the neutrino 
energy deposition behind the shock. 

While the shock is stalled, matter continues to accrete on the proto-neutron star 
passing through the shock. During this so-called “accretion phase” the shocked hot 
material behind the shock, composed mainly by free nucleons, electrons and photons, is 
heated by the accretion and produces neutrinos and antineutrinos of all flavors through 
the processes 4H3MHD. Since the stalled shock is out of the neutrinosphere, these 
neutrinos can free-stream out of the star and cause the so-called “hump” in the neutrino 

1( The recoil kinetic energy of the final proton is negligible. Indeed, momentum conservation implies 
that the momentum p p of the final proton is of the order of the momentum p Ue of the initial neutrino, 
which is practically equal to the neutrino energy, because of the smallness of neutrino masses. Since the 
neutrino energy is smaller than a few tens of MeV, the recoil kinetic energy of the proton, Pp/2m p , is 
suppressed by the large mass rri p of the proton. 
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Figure 17: Time evolution of neutrino luminosity and average energy obtained with the 
numerical supernova model in Ref. |2b6] . The dashed lines are for zz e , the solid lines for 
P e , and the dot-dashed lines for i/ T , P r , which are collectively denoted by v x . The 

neutronization burst is visible as a peak of luminosity and energy of electron neutrinos, 
which happens 3-4 msec after bounce. 


luminosity curve shown in Fig. El The average neutrino energy is low during the hump 
because the dense matter in the shock is opaque to high-energy neutrinos. As the shock 
gradually revives, the matter density decreases and the average neutrino energy increases. 

Summarizing, in the prompt explosion scenario there are two phases of the neutrino 
flux: first a brief and intense burst of prompt electron neutrinos from shock breakout, 
with a degenerate spectrum of high energy, which is however so brief that little energy 
(about 10 51 ergs) and lepton number are carried away. Then there is a less intense thermal 
emission of neutrinos of all flavors which last for a few seconds and carries away most of 
the binding energy of the neutron star (about 3 x 10 53 ergs). The total number of emitted 
neutrinos and antineutrinos exceeds by an order of magnitude the original lepton number 
of the collapsed core. 

In the delayed explosion scenario, in addition to the prompt electron neutrino burst 
and the thermal emission of neutrinos of all flavors one expects an accretion phase which 
prolongs the peak of the thermal neutrino luminosity over a time scale of about half 
second. 

The delayed explosion scenario constitutes a sort of standard model of core-collapse 
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supernova explosion. However, the possibility of shock revival through neutrino heating 
is still under study (see Ref. fznaESHi). 

8.3 SN1987A 

On 24 February 1987 a new very bright Type II supernova, SN1987A, was discovered in 
the Large Magellanic Cloud, which is a satellite galaxy of the Milky Way, at a distance 
of about 50 kpc from the solar system (see Refs. [K31I54 ] ). At that time four large under¬ 
ground neutrino detectors potentially sensitive to supernova neutrinos were in operation: 
Kamiokande-II [26411265] . IMB [266U267U268] . Baksan [2691127011271] and LSD |272j . These 
detectors saw an unusual number of events with energy of the order of 10 MeV within a 
time window of the order of 10 sec in the hours before the optical discovery of SN1987A. 
The events observed in the Kamiokande-II, IMB and Baksan happened at the same time 
(within uncertainties of the absolute time calibration of the detectors and the random 
occurrence of the events), whereas the LSD events have been recorded about five hours 
before those of the other detectors, at a time when the other detectors did not see any sig¬ 
nal. Therefore, there is a controversy on the origin of the LSD events (see Refs. |2ZnU2ZU) 
and usually the LSD events are not included in the analysis of SN1987A data. In the 
following sub-subsections we describe briefly the data of the Kamiokande-II, IMB and 
Baksan detectors, which are used to set limits on neutrino properties. 

Supernova SN1987A is the best studied of all supernovae not only because of the 
detection of its neutrinos but also because it was the first supernova visible to the naked 
eye after the Kepler in 1604 and because it is the only supernova for which the progenitor 
star is known: it was a blue supergiant B3 I star named Sanduleak - 69°202 |ZZS|. 

The evolution of the remnant of SN1987A has been deeply studied in all spectral 
bands (see references in Refs. EH El ED- Although no compact remnant has been 
identified with certainty so far, there is some indication of the presence of a 2.14 ms 
optical pulsar [27iB j. 

The observation of SN1987A neutrinos marked the beginning of extra solar system 
neutrino astronomy 18 . It has been one of the great achievements of the Kamiokande 
detector, which was designed by Masatoshi Koshiba and earned him the 2002 Nobel 
Prize in Physics. 

8.3.1 Kamiokande-II 

The Kamiokande detector (see [2771I27E1 ) was a water Cherenkov detector with a fiducial 
volume containing 2140 tons of water surrounded by 948 photomultiplier tubes of 50 cm 
diameter, covering about 20% of the surface area. It was located in the Kamioka mine in 
Gifu prefecture, Japan, with a 2400 rn.w.e. overhead shielding. 

The Kamiokande detector was built in 1983 for the search of nucleon decay (Kamiokande 
is the acronym of Kamioka Nucleon Decay Experiment), although the possible detection 
of supernova neutrinos was mentioned in the original proposal (see Ref. [277] )- In 1986 the 
Kamiokande detector was upgraded to Kamiokande-II for the detection of solar 8 B neutri¬ 
nos with a threshold of about 6 MeV. In 1990 the detector was upgraded to Kamiokande- 

18 Solar neutrino astronomy was started by Raymond Davis Jr., co-winner of the 2002 Nobel Prize in 
Physics, with the Homestake Chlorine experiment. 


62 







Ill and continued operation until 1995. Besides the search for nucleon decay P79ll280j and 
the observation of solar 8 B neutrinos BED, the Kamiokande detector obtained two unex¬ 
pected important results during the Kamiokande-II phase: the observation of SN1987A 
neutrinos [264112651] and the discovery of the atmospheric neutrino anomaly [I282li283| . In 
1996 the Kamiokande detector was replaced by the Super-Kamiokande detector located 
in the same mine, which has a fiducial volume of 22.5 kt (see Ref. [ 284 ]). 

After the optical discovery of supernova SN1987A the Kamiokande-II collaboration 
examined carefully their data looking for a significant number of events above background 
in a time interval of the order of 10 sec and energy of the order of 10 MeV. They found 
such a collection of events at 7:35:35 UT of 23 February 1987. Unfortunately, before the 
discovery of supernova SN1987A the Kamiokande-II collaboration did not think that an 
accurate time measurement was necessary and the clock of the experiment was set by 
hand. As explained in Ref. [ 265] . “it would be straightforward after SN1987A to have 
made an absolute calibration of the clock ..., but and abrupt power outage took place 
in the Kamioka mine on 25 February 1987, and precluded that alternative measure”. 
Therefore there is an uncertainty of about one minute in the Kamiokande-II determination 
of the time in which the SN1987A neutrino burst passed the Earth. 

Electron antineutrinos with energy larger than 1.8 MeV can be detected with the 
“inverse /3-decay” reaction 

P e + p —> n + e + . (8.13) 


with cross section 

f E \ 2 

a De+p ^n+e+ - 8.5 x 1(T 44 cm 2 . (8.14) 

The produced positron is emitted almost isotropically and can be observed in water 
Cherenkov detectors, as well as in scintillator detectors. The energy of the incident v e is 
given by 

E Pe = E e + + m n — m p = E e + + 1.3 MeV . (8.15) 

The Kamiokande-II detector observed also neutrinos through the elastic scattering 
reaction 

v T e — —> v T e~ , (8.16) 


which is mainly sensitive to electron neutrinos (indeed, it has been used to observe solar 
electron neutrinos): 


a 


ES 


~3<t£~6<t5, 


(8.17) 


with 

~9.2xl(T 45 2Kcm 2 . (8.18) 

Table El shows the relative time t, the energy E e and the angle 9 e of the observed 
charged lepton with respect to the direction opposite to SN1987A of 12 events measured 
in the Kamiokande-II detector during the supernova SN1987A neutrino burst. The event 
number 6 is reported in the Kamiokande II original publication 1265] . but is excluded in 
their signal analysis because its low number of hit photomultipliers indicates that it is 
likely to be due to background. 

Most authors agree that it is most likely that all Kamiokande events have been gen¬ 
erated through the inverse /3-decay reaction dHU because of the dominance of 
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Kamiokande II 


Event 

Time t 
(sec) 

Energy E e 
(MeV) 

Angle 9 e 
(degrees) 

1 

0 

20.0 ±2.9 

18 ± 18 

2 

0.107 

13.5 ±3.2 

40 ±27 

3 

0.303 

7.5 ±2.0 

108 ±32 

4 

0.324 

9.2 ±2.7 

70 ±30 

5 

0.507 

12.8 ±2.9 

135 ± 23 

6 

0.686 

6.3 ± 1.7 

68 ±77 

7 

1.541 

35.4 ±8.0 

32 ± 16 

8 

1.728 

21.0 ±4.2 

30 ± 18 

9 

1.915 

19.8 ±3.2 

38 ±22 

10 

9.219 

8.6 ±2.7 

122 ± 30 

11 

10.433 

13.0 ±2.6 

49 ±26 

12 

12.439 

8.9 ± 1.9 

91 ±39 


Tabic 6: Relative time t, energy E e and angle 6 e of the observed charged lepton with 
respect to the direction opposite to SN1987A of the Kamiokande II events [2ESj ■ The 
event number 6 is likely to be due to background because it has a low number of hit 
photomultipliers. 


its cross section. Nevertheless, some authors 265. 286; have speculated on the fact that 
the first event point almost in the opposite direction of the LMC 19 , which could be an 
indication that it is due to an electron neutrino interacting in the detector trough the 
elastic scattering process dump - 

8.3.2 IMB 

The IMB water Cherenkov detector was located in a salt mine near Fairport, Ohio, USA, 
at a depth of 1570 rn.w.e.. It consisted of a rectangular tank filled with purified water with 
an active volume of about 6800 tons viewed by 2048 8-inch photomultipliers arranged on 
an approximate 1 m grid. 

On 23 February 1987 the IMB detector recorded eight neutrino-produced events with 
energies between 20 and 40 MeV in a time interval of 6 sec starting from 7:35:41.37 UT 
(the clock had an absolute uncertainty of 50 msec and a relative uncertainty of 0.5 msec). 
The background rate is negligible, about 2 per day in the range 20-2000 MeV. 

The important characteristics of the eight IMB events are listed in Tabled Since these 
events are most likely due to the inverse /3-decay process turn the neutrino energy is 
given by Eq. lETTSt . 

Taking into account the trigger efficiency and about 13% dead time of the detector, 
the IMB collaboration estimated that 35 ± 15 neutrino events with energy above 20 MeV 
occurred in the detector | 22 !. 

19 In the first Kamiokande II publication, Ref. |2d4j , the angle of the second event was reported to be 
15 ± 27, pointing almost backward from the direction of the Large Magellanic Cloud. This angle was 
corrected in the second Kamiokande II publication, Ref. EM]. 
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IMB 

Event 

Time t 
(sec) 

Energy E e 
(MeV) 

Angle 9 e 
(degrees) 

1 

0 

38 ±7 

80 ± 10 

2 

0.412 

37 ±7 

44 ± 15 

3 

0.650 

28 ±6 

56 ±20 

4 

1.141 

39 ±7 

65 ±20 

5 

1.562 

36 ±9 

33 ± 15 

6 

2.684 

36 ±6 

52 ± 10 

7 

5.010 

19 ±5 

42 ±20 

8 

5.582 

22 ±5 

104 ± 20 


Tabic 7: IMB supernova SN1987A events from Ref. (26 7j . The time of each event is 
relative to the first one, which occurred at 7:35:41.37 UT of 23 February 1987. There is 
an additional systematic uncertainty in the energy scale estimated to be about 10%. The 
background rate is negligible (about 2 events per day). 

8.3.3 Baksan 

The Baksan Underground Scintillation Telescope 209 . 270 271 is located in the Baksan 
neutrino Observatory at a depth of 850 m.w.e. in the Baksan Valley in North Caucasus, 
Russia. The telescope consists of 3150 parallclepipedal tanks filled with oil-based liquid 
scintillator viewed by a 15 cm photomultiplier. The energy threshold for supernova 
neutrinos is about 10 MeV. The total target mass is about 330 tons. The background, 
mainly caused by cosmic ray muons and discharges in the photomultipliers, is relatively 
large. Therefore, only 1200 inner tanks with lower background and a mass of about 130 
tons are used as signal triggers, and the inner tanks plus part of the external tanks are 
used as fiducial volume, with a mass of about 200 tons. 

As water Cherenkov detectors, the Baksan Underground Scintillation Telescope is 
mostly sensitive to electron antineutrinos which interact with protons through the inverse 
/1-decay reaction dHUSD - 

At the time of SN1987A the Baksan Underground Scintillation Telescope had been in 
operation for about six years. During this period of time, including 23 February 1987, it 
never happened that more than 7 events were observed in an interval of 20 sec. The Bak¬ 
san Collaboration were expecting about 35 antineutrino events in the trigger mass and 
about 54 events in the fiducial mass for a supernova at a distance of 10 kpc (i.e. within 
the Milky Way). In the period from 1 to 23 February 1987 the Baksan Underground 
Scintillation Telescope did not measure pulse clusters that differ significantly from the 
background. Therefore, the Baksan Collaboration could not claim an independent ob¬ 
servation of SN1987A neutrinos. However, when supplemented by the information of the 
Kamiokande-11 and IMB observations, the Baksan Collaboration identified 5 events in a 
10 sec interval that may overlap with the Kamiokande-II and IMB, taking into account an 
uncertainty of 7g 4 sec in the absolute Baksan clock measurement. The Baksan clock had 
a relative accuracy of about one millisecond and a nominal absolute accuracy of about 
2 sec, but on 11 March 1987 it was found that the clock had developed a forward shift 
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Baksan 

Event 

Time t 

(sec) 

Energy E e 
(MeV) 

1 

0 

12 ±2.4 

2 

0.435 

18 ±3.6 

3 

1.710 

23.3 ±4.7 

4 

7.687 

17 ±3.4 

5 

9.099 

20.1 ±4.0 


Tabic 8: Relative time t and energy E e of the Baksan SN1987A events 1221- 

in time of 54 sec that conld have happened in one step or gradually since 17 February 
1987. Since the Baksan clock time of the five Baksan events is about 30 sec after the 
IMB events (which were measured with absolute time uncertainty of about 50 msec), the 
simultaneous occurrence of Baksan and IMB events is possible. 

Table |H1 shows the relative time t and the energy E e of the five Baksan events |27nj . 
However, since the background rate in the Baksan detector is rather high, it is impossible 
to know which, if any, of the events is due to SN1987A neutrinos. Therefore, most authors 
did not include the Baksan data in the analysis of SN1987A neutrino events. 

8.3.4 Comparison with Theory 

The neutrino events have been compared with theoretical predictions in many papers 
|287U288ll285LI289ll274U2901 l 291LI252| . Although only about two dozens of the estimated 
10 28 neutrinos that passed through the Earth were detected, these few events delivered us 
precious information about the physics of core-collapse supernovae. Most authors agree 
that the detected neutrino events are compatible with the general features of the standard 
core-collapse supernova scenario described in section IH21 

The most accurate analysis of SN1987A neutrino data has been performed recently 
by Loredo and Lamb |252j , which, for the first time, took into account the background in 
the Kamiokande-II and Baksan detectors. This is important, because it is impossible to 
know with certainty which events have been really produced by neutrinos coming from 
SN1987A and which events are due to background. 

Table El shows the relative time t, the energy E e , the estimated background rate 
B(E e ), and the probabilities Pb (prompt) and Pb (delayed) that the event is due to back¬ 
ground according to the best fit prompt and delayed supernova explosion models (see 
section I8.2j) of the 16 Kamiokande-II events taken into account in the analysis of Loredo 
and Lamb !2E2J- The events number 13-16 have not been considered as SN1987A events 
by the Kamiokande II Collaboration (see Table El), although they can be seen in Fig. 9 of 
Ref. . Indeed, from the last two columns of Table El one can see that, according to the 
calculation in Ref. |252j . these events have a high probability to be due to background. 
However, the probability that at least one of them is a signal event is not negligible and 
it is correct to include them in the data analysis, as done in Ref. [222!- 

From Table El one can also see that the event number 6, which was excluded from 
the Kamiokande-II signal analysis [ 265| has indeed a non-negligible probability to be 
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Kamiokande 11 

Event 

Time t 
(sec) 

Energy E e 
(MeV) 

B(E e ) (252; 
(s- 1 ) 

Pb (252^ 
(prompt) 

Pb j252] 
(delayed) 

1 

0 

20.0 ±2.9 

1.6 x 10~ 5 

5.8 x 10" 5 

2.4 x 10~ 5 

2 

0.107 

13.5 ±3.2 

1.9 x 10~ 3 

6.3 x 10" 3 

1.9 x 10~ 3 

3 

0.303 

7.5 ±2.0 

2.9 x 10~ 2 

0.16 

4.7 x 10- 2 

4 

0.324 

9.2 ±2.7 

1.2 x 10" 2 

5.4 x 10- 2 

1.7 x 10- 2 

5 

0.507 

12.8 ±2.9 

2.1 x 10~ 3 

7.6 x 10" 3 

3.2 x 10~ 3 

6 

0.686 

6.3 ±1.7 

3.7 x 10~ 2 

0.25 

0.15 

7 

1.541 

35.4 ±8.0 

4.5 x 10~ 5 

1.2 x 10“ 3 

1.5 x 10~ 3 

8 

1.728 

21.0 ±4.2 

8.2 x 10~ 5 

5.7 x 10" 4 

1.0 x 10~ 3 

9 

1.915 

19.8 ±3.2 

1.5 x 10~ 5 

9.9 x 10" 5 

1.9 x 10~ 4 

10 

9.219 

8.6 ±2.7 

1.5 x 10~ 2 

0.33 

0.49 

11 

10.433 

13.0 ±2.6 

1.9 x 10“ 3 

0.11 

0.12 

12 

12.439 

8.9 ± 1.9 

1.6 x 10~ 2 

0.54 

0.60 

13 

17.641 

6.5 ± 1.6 

3.8 x 10~ 2 

0.92 

0.89 

14 

20.257 

5.4 ± 1.4 

2.9 x 10~ 2 

0.97 

0.94 

15 

21.355 

4.6 ± 1.3 

2.8 x 10" 2 

0.97 

0.93 

16 

23.814 

6.5 ± 1.6 

3.8 x 10~ 2 

0.99 

0.94 


Tabic 9: Relative time t, energy E e , event background rate B(E e ), and probabilities 
Pb (prompt) and Pg (delayed) that each event is due to background in the best fit prompt 
and delayed supernova explosion models of the Kamiokande-II events taken into account 
in Ref. 252j . 


a background event according to the best fit prompt and delayed supernova explosion 
models (see section 18.21) calculated in Ref. 12521 . 

The ability to take into account background events of the Loredo and Lamb method 
[2521 is mostly useful for the inclusion in the analysis of SN1987A of the Baksan data. 
Table imi shows the relative time t, the energy E e , the event background rate B(E e ), and 
probabilities Pb (prompt) and Pb (delayed) that each event is due to background in the 
best fit prompt and delayed supernova explosion models in Ref. 252 of the Baksan events. 
One can see that the background rate in the Baksan detector is rather high. For this 
reason most authors did not include the Baksan data in the analysis of SN1987A neutrino 
events. However, Loredo and Lamb ^2521 properly took into account the background 
rate and proved that the Baksan events are compatible with a supernova signal. The 
probabilities Pb (prompt) and Pb (delayed) show that some of the Baksan events could 
be due to supernova electron antineutrinos. 

Loredo and Lamb ( 2521 found that models of supernova explosion with the delayed 
mechanism explained in section 18.21 are about 100 times more probable than prompt 
explosion models. The electron antineutrino average energy is about 15 MeV, as expected 
from the cooling of the proto-neutron star (see Eq. (JHD2D). The cooling time scale is about 
4 sec, and the time scale of the accretion component is about 0.7 sec, in agreement with 
numerical calculations. The total inferred number of electron antineutrinos emitted is 
about 3 x 10 57 , implying a binding energy of the neutron star of about 3 x 10 53 ergs, 
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Baksan 

Event 

Time t 

(sec) 

Energy E e 
(MeV) 

/i (/•;,) |2!2] 
(s- 1 ) 

Pb |252| 

(prompt) 

Pb f252] 
(delayed) 

1 

0 

12.0 ±2.4 

8.4 x 10~ 4 

2.1 x 10“ 2 

4.9 x 10~ 3 

2 

0.435 

17.9 ±3.6 

1.3 x 10~ 3 

3.6 x 10- 2 

1.9 x 10- 2 

3 

1.710 

23.5 ±4.7 

1.2 x 10~ 3 

7.4 x 10" 2 

0.12 

4 

7.687 

17.6 ±3.5 

1.3 x 10~ 3 

0.30 

0.35 

5 

9.099 

20.3 ±4.1 

1.3 x 10~ 3 

0.55 

0.52 


Table 10: Relative time t, energy E e , event background rate B(E e ) and probabilities 
Pb (prompt) and Pb (delayed) that each event is due to background in the best fit prompt 
and delayed supernova explosion models of the Baksan events taken into account in 
Ref. j2S2| . 


as expected from simple estimations (see section E2D- Unfortunately, as explained in 
Ref. j.22j, the SN1987A neutrino data are too sparse to obtain more detailed information 
on the supernova mechanism. 


8.4 Neutrino Mass 

The basic idea of constraining neutrino masses from the observation of supernova neutri¬ 
nos was proposed many years ago |551l5ifllI)?U58ll5§] . 

An extremely relativistic neutrino with mass m E propagates with a group velocity 



(8.19) 


If a neutrino flux is emitted by a source at a distance D, the time-of-flight delay of a 
massive neutrino with respect to a massless particle (as a photon or a graviton) emitted 
by the same source is 


At 



v 


771 

2TP 


D 



E V 2 D 

- 1 -sec 

MeV ) 50kpc 


( 8 . 20 ) 


If neutrinos are emitted in a burst with intrinsic duration AT 0 , the observation of events 
separated by a time interval larger than ATq would provide a direct measurement of the 
neutrino mass (assuming D known and E measurable). If the neutrino energy spectrum 
has mean value E and width AE, neutrinos produced at the same time with different 
energies would reach a detector at a distance D in the time interval 


AT ~ 


m 2 AE 
EU “ft" 


( 8 . 21 ) 


The model-independent sensitivity to the neutrino mass is found by requiring this time 
interval to be smaller than the intrinsic duration of the neutrino burst: 


AT < AT 0 < AT obs , 


( 8 . 22 ) 
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Authors 

Upper Bound 

Bahcall and Glashow (1987) [292] 
Arnett and Rosner (1987) 293 

Sato and Suzuki (1987) |285j 
Schramm (1987) [274] 

Krauss (1987) (2T?Hj 

Spergel and Bahcall (1988) [2511 
Abbott, De Rujula and Walker (1988) [295j 
Loredo and Lamb (2001) .252) 

?% e < 11 eV 
m Pe < 12 eV 
m Pe < 26 eV 
m Pe 7: 30 eV 
m Pe < 8 - 15 eV 
m Pe < 16 eV (95% C.L.) 
rn Pe < 11 eV (95% C.L.) 
m 9e < 5.7eV (95% P.) 


Table 11: Upper bounds for m Pe obtained by several authors from the analysis of 
SN1987A neutrino data. C.L. = Confidence Level in Frequentist analyses and P. = 
Probability in Bayesian analyses. 


where AT 0 b s is the observed time interval of the neutrino burst. The inequalities (18.221) 
imply the upper bound 


m < E 


E A T t 


obs 


A E D 


~ 14 eV 


10 MeV 


E ( AT C 


obs 


A E \ 10 sec 


1/2 


50 kpc 
D 


1/2 


(8.23) 


It is clear that a large distance, a quick neutrino burst, a low neutrino energy, and a 
wide energy range are advantageous for the measurement of an effect due to the neutrino 
mass. Unfortunately, increasing the distance decreases the neutrino flux at the detector 
in proportion to D~ 2 and decreasing the energy decreases the detection event rate. In 
practice, the energy of neutrinos coming from a supernova is of the order of 10 MeV 
and the existing detectors allow to observe only neutrinos produced by supernovae in our 
galaxy or in its satellites (the Small and Large Magellanic Clouds). 

Supernova SN1987A occurred in the Large Magellanic Cloud, at a distance of about 50 
kpc from the Solar System. Since the measured neutrino burst had an average energy E ~ 
15 MeV, a width A E ~ 15 MeV, and an estimated original time duration AT 0 b s ~ 10 sec, 
from Eq. ([8.23]) one can see that the observation of the neutrinos SN1987A allows a model- 
independent sensitivity to a neutrino mass m > 20 eV. Since the time duration of the 
neutrino signal is compatible with theoretical predictions, the SN1987A data allow only to 
obtain an upper limit on the neutrino mass of the order of 20 eV. Indeed, Schramm mu 
argued that without making specific model assumptions, all that can be safely said is 


rri L 


<30eV. 


(8.24) 


Many authors have calculated an upper bound on the electron antineutrino mass 
from the SN1987A neutrino data with some specific assumptions, often well-motivated, 
about the intrinsic spread of the neutrino burst, obtaining upper bounds for m Pe lying 
in the 5 — 30 eV range |292[l293[l285[l274M290M294ll295M252| , as shown in Table HU These 
bounds were also obtained with different statistical techniques for the analysis of the few 
available events. However, Loredo and Lamb [ 252) noticed that not all of these statistical 
procedures are appropriate. 

In their accurate recent analysis of SN1987A neutrino data, Loredo and Lamb mu 
applied the Bayesian method, which is rather easily implemented in a correct way and 
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leads to results with a clear meaning (see also Ref. 2961 1. Their upper limit on the 
electron antineutrino mass is 

m Pe < 5.7eV, (8.25) 

with 95% probability. This limit is comparable with the other most stringent existing 
limits on the electron neutrino mass. It is significantly more stringent than Schramm’s 
model-independent limit (I8.29jl because it follows from a detailed fit of the data in terms 
of a delayed explosion model, which is favored by the data as explained in section 18.5.41 

8.4.1 Neutrino Mixing 

If there is neutrino mixing (see Section 0 ), an electron antineutrino does not have a 
definite mass, since it is a superposition of different massive neutrinos. In this case the 
kinematical upper limit (18.2511 applies to all the massive neutrinos that have a substantial 
mixing with D e . 

Current experimental data on solar and atmospheric neutrinos indicate the existence 
of three-neutrino mixing (see Section 0 with 

Am 2 ol « 5 x 10 " 5 eV 2 , Am 2 tm « 2.5 x 10 ~ 3 eV 2 , (8.26) 

where Am^ and Am 2 tm are, respectively, the effective neutrino squared-mass differences 
in two-neutrino analyses of solar and atmospheric data. Since the squared-mass differ¬ 
ences in Eq. (18.2611 are very small, the kinematical upper limit (18.2511 applies to all the 
three neutrino masses mi, m 2 and m 3 . 

Since the Kamiokande-II SN1987A events appear to be clustered in time in two groups 
separated by an interval of about 10 sec, some authors [237112981 have claimed that there 
is an evidence of two mass groupings at about 4 eV and 22 eV 2TJ7. However, these 
authors had to assume that electron antineutrinos are emitted from the supernova in a 
very short time, of the order of 0.1 sec. This assumption is contrary to our understanding 
of the core-collapse supernova mechanism, according to which electron antineutrinos are 
emitted during the cooling phase of the proto-neutron star on a time scale of about 10 
sec (see section 18.2(1 . Moreover, the existence of neutrinos with masses of about 4 eV and 
22 eV which have large mixing with the electron antineutrino is excluded by the Tritium 
upper bound on the effective electron antineutrino mass (see Section 0 . 

Other information on neutrino mixing has been obtained from SN1987A data consid¬ 
ering the effect of vacuum oscillations or MSW mm resonant transitions on the fluxes 
of different flavors. Large v x v e transitions are disfavored, because they would imply 
a harder spectrum of h e ’s on Earth than observed (see Refs. [2991130111180 lj and references 
therein). 

8.4.2 Other Neutrino Properties 

For the sake of completeness, let us briefly list some of the other neutrino properties that 
have been constrained using SN1987A neutrino data. 

Since electron antineutrinos arrived at the Earth from a distance of about 50 kpc, 
their lifetime r Pe is constrained by [2EK1I2731J 

Re £ 1-6 X 10 5 {m 0e /E v J yr. (8.27) 
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The total amount of emitted energy inferred from the measured v e flux is compatible 
with the binding energy of a neutron star only if the number N v of neutrino flavors is 
limited by |3021l274l f290j 

K<6. (8.28) 

The cooling of the proto-neutron star constrains the Dirac masses of and u T by 
}3n3ii3n4il3n5llMl3n7j 


m“ rac < 14keV 

Vfj, 

m Di rac < 14 keV or m°; rac >34MeV. (8.29) 

The absence of 7 emission accompanying the SN1987A neutrino burst implies a lower 
bound between about 10 6 and 10 10 yr for the lifetime of a heavy massive neutrino with 
mass 2 m e < rrih < 100 MeV which as a substantial mixing with the active light flavor 
neutrinos and decays via u h —■► v k + e + + e~ |308U309U310j . 

The observed 10 sec timescale of cooling of the proto-neutron star implies an upper 
bound mEsmmm- 

^<10~ 12 ^ b (8.30) 

for the electron neutrino magnetic moment which could flip neutrino hclicity through scat¬ 
tering with electrons and nucleons or through interactions with the strong magnetic field, 
generating sterile right-handed neutrinos that escape freely, cooling the proto-neutron 
star in less than 1 sec. 

The absence of a similar cooling by right-handed neutrino emission constrains also 
the charge radius of right-handed neutrinos by M 

(r 2 ) R < 2 x 10' 33 cm 2 . (8.31) 

The electric charge of the electron neutrino is bounded by EE! 

<!,, < 1" 'V (8.32) 

otherwise the galactic magnetic field would have lengthened the neutrino path and neu¬ 
trinos of different energy could not have arrived on Earth within a few seconds. 

8.5 Future 

Several detectors sensitive to supernova neutrinos are currently in operation (Super- 
Kamiokande |2E3|, SNO HHZj, LVD [HU ], KamLAND 0, AMANDA [SH] , MiniBooNE 
|320| ) or under preparation or study (see Ref. [42] and references therein). Many au¬ 
thors have studied future possibilities of supernova neutrino detection and its potential 
sensitivity to neutrino masses (see Refs. [250.. 1321, I5T2 . f.T2Tl .j and references 
therein). Current and future supernova neutrino detectors are much larger than the de¬ 
tectors in operation during 1987 and the order of magnitude of the total number of events 
expected when the next galactic supernova will explode is 10 J . Such impressive statistics 
will be precious in order to test our understanding of supernova physics and improve our 
knowledge of neutrino properties. 
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There is a general agreement among workers in the field that future supernova neutrino 
detections cannot be sensitive to an effective 20 electron neutrino mass smaller than a few 
eV, because of the intrinsic spread in time of the neutrino burst. Totani 222 has shown 
that using the correlation between neutrino energy and arrival time implied by Eq. flQUP , 
it is possible to reach a sensitivity of about 3 eV for the effective electron neutrino mass. 
Beacom, Boyd and Mezzacappa [325113291 have shown that an abrupt termination of the 
neutrino signal due to black-hole formation may allow the Super-Kamiokande detector to 
be sensitive to an electron neutrino mass as low as 1.8 eV. Another interesting possibility 
is the measurement of the time delay between gravitational waves generated by core 
collapse and the neutronization neutrino burst ]531l591l327j , which may allow a sensitivity 
to the neutrino mass of about 1 eV [32%). 

However, since the current upper limit for the effective electron neutrino mass is 
already a few eV (see Section 01) and the future KATRIN experiment [2B.] will be able 
to push the limit down to about 0.3 eV, a supernova limit on m Ve will not be extremely 
exciting. Therefore, several authors have concentrated on the possibility to constrain the 
effective masses of and u T j25(HI321ll323ll324ll32f)ll32d] . whose laboratory limits are well 
above the eV scale (see Section 01). 

The flux of supernova z^, v T and P T is of the same order as that of v e and P e , but 
the problem is to distinguish them, because they can be observed only through neutral- 
current interactions, which are flavor blind (the energy is too low to produce p or r in 
charged-current reactions). Therefore, the u T , v r signal can be only extracted on 

a statistical basis by subtracting the v e and D e contributions from the measured neutral- 
current signal. The v e and v e contributions are estimated from the v e and v e charged- 
current signals. 

Unfortunately, in usual neutral-current neutrino interactions, as that in SNO, 

v + d —>p + n + z/, (8.33) 

the energy of the neutrino is not determined. Therefore, it is not possible to use the 
correlation between neutrino energy and arrival time implied by Eq. flHUDt for the mea¬ 
surement of neutrino masses, and the upper limit on the effective masses of and u T 
cannot be pushed below about 30 eV mmm- An interesting exception is the 
abrupt termination of the neutrino signal due to black-hole formation, which may allow 
a sensitivity to and u T effective masses as low as about 6 eV Another 

promising technique which has been recently proposed by Beacom, Farr and Vogel [329 ] 
is the measurement of the recoil proton kinetic energy in neutral-current neutrino-proton 
elastic scattering, 

v + p —> v + p . (8.34) 

Since the recoil protons have a kinetic energy of the order of 1 MeV, they are non- 
relativistic and cannot be seen in water Cherenkov detectors, but they can be observed 
in liquid scintillator detectors as KamLAND j9) and BOREXINO [330] . Unfortunately, 
the proton direction cannot be measured in scintillator detectors, denying the possibility 
to reconstruct the neutrino energy from simple kinematics on a event-by-event basis. 

20 In this context we use the adjective “effective” for the masses of flavor neutrinos in order to keep 
in mind that in reality what are measured are the masses of the massive neutrinos which have a large 
mixing with the flavor neutrino under consideration. 
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However, Beacom, Farr and Vogel ji3'29) have shown that a fit of the proton kinetic energy 
distribution could allow to measure the neutrino temperature and the total neutrino 
energy with an accuracy of about 10%. As far as we know, the possibility to obtain 
information on the effective masses of and v T with this method has not been explored 
so far. 

Of course a major problem in supernova neutrino physics is the actual occurrence of 
a supernova at a galactic scale distance. As we have seen in section 18.11 the estimated 
rate of core-collapse supernovae in the Milky Way is about 2 ± 1 per century. Such low 
rate is just at the border of the patience of very patient scientists. Since most scientists 
are not so patient, there is an active research to study the feasibility of huge detectors 
that could observe a few tents of events produced by a supernova in the local group of 
galaxies (see Ref. and references therein). 

9 Conclusions 

During many years there were indications in favor of the disappearance of solar v e 's (the 
so-called “solar neutrino problem”), obtained in the Homestake H2|. Kamiokande (2SU, 
GALLEX [O] and SAGE [HZ] solar neutrino experiments, and indications in favor of 
the disappearance of atmospheric z/^’s (the so-called “atmospheric neutrino anomaly”), 
found in the Kamiokande |331M332ll283| . IMB |333j . Soudan 2 [10J and MACRO [TT] 
atmospheric neutrino experiments. Neutrino oscillations were considered to be the natural 
interpretation of these data, in spite of the fact that other possibilities were not excluded 
at that time (as large anomalous magnetic moment of neutrino, neutrino decay, etc.; see 

Refs. |3341l335llMi251 ). 

In the latest years with the impressive results of the atmospheric and solar neutrino 
experiments Super-Kamiokande and SNO and the long-baseline reactor experiment Kam- 
LAND, the status of neutrino masses, mixing and oscillations drastically changed. The 
up-down asymmetry of the atmospheric multi-GeV muon events discovered in the Super- 
Kamiokande experiment jT] , the evidence of transitions of solar v e into v M and v T obtained 
in the SNO experiment [7] from the observation of solar neutrinos through the detection 
of CC and NC reactions, and the evidence of disappearance of reactor v e found in the 
long-baseline KamLAND experiment |2] imply that neutrino oscillations driven by small 
neutrino masses and neutrino mixing is the only viable explanation of the experimental 
data. 

The generation of neutrino masses, many orders of magnitude smaller than the masses 
of their family partner leptons and quarks, requires a new mechanism beyond the Stan¬ 
dard Model. Several possibilities for new mechanisms which could generate small neutrino 
masses are open today. They are based on the see-saw type violation of the total lepton 
number at a very large scale, large extra dimensions, etc. (see Refs. B3SI3321I33HIIBH). 

It is obvious that the understanding of the true mechanism of the generation of neutrino 
masses and mixing requires new experimental data. 

First of all, we need to know how many massive neutrino exist in nature. The mini¬ 
mal number of massive neutrinos is equal to the number of flavor neutrinos (three). If, 
however, the LSND indication in favor of short-baseline > v e oscillations is confirmed, 
at least four massive neutrinos are needed. In spite of the fact that such scenario is 
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disfavored by the analysis of the existing data, in order to reach a conclusion it is needed 
to wait the decisive check of the LSND result which is under way at Fermilab with the 
MiniBooNE experiment. 

The problem of the nature of the massive neutrinos (Dirac or Majorana?) is crucial 
for the understanding of the origin of neutrino masses. This problem can be solved by the 
experiments searching for neutrinoless double-/? decay of some even-even nuclei. From the 
existing data for the effective Majorana mass \(m) \ the bound \(m)\ < (0.3 — 1.2) eV has 
been obtained. Several new experiments with a sensitivity to \(m)\ improved by about 
one order of magnitude with respect to present experiments are under preparation. 

Neutrino oscillation experiments allow to determine the neutrino mass-squared differ¬ 
ences Am 2 . Since neutrino oscillations is an interference phenomenon, neutrino oscillation 
experiments are sensitive to very small values of Am 2 . On the other hand, the problem 
of the absolute values of neutrino masses is apparently the most difficult one in the ex¬ 
perimental investigation of the physics of massive and mixed neutrinos. In this paper we 
presented a review of our present knowledge of the absolute values of neutrino masses 
and the prospects for the future (see also Ref. (22). 

At present, the bound mp < 2.2 eV at 95% C.L. has been obtained from /?-decay 
experiments jHUJESl- About one order of the magnitude improvement is expected in the 
future. 

With enormous progress in cosmology in the last years the possibility to obtain infor¬ 
mation about the values of neutrino masses from cosmological data has strongly increased. 
The existing data allowed to obtain the bound rn t < 2 eV for the total mass of neu¬ 
trinos mmm- About one order of magnitude improvement is expected with the 
future data. 

A possible explanation of the observation 

of very high-energy cosmic rays beyond the GZK cutoff mm is the existence of Z- 
bursts which are possible if there is at least one neutrino mass in the 0.08 — 1.3 eV 

range at 68% C.L.. 

The observation of supernova 1987A neutrinos by the Kamiokande-II [265j. IMB nm 
and Baksan [27T)| detectors allowed to obtain important information on supernova physics 
and neutrino properties. The absence of an energy-dependent spread of the neutrino burst 
constrains the absolute value of the effective electron neutrino mass below 5.7 eV with 
95% probability 2o2 . 

In the framework of the minimal scheme with three massive neutrinos the values of 
the neutrino masses are determined by the minimal neutrino mass mi and the neutrino 
mass-squared differences Aand Arn 2 2 , which will be measured with high precision 
in solar, atmospheric, long-baseline reactor and accelerator neutrino oscillations experi¬ 
ments. Thus, the problem of the absolute values of neutrino masses is the problem of 
determination of the minimal mass m\. If m\ is much smaller than Am 2 ol ~ 7 • 10 -3 eV 
(hierarchical mass spectrum), it will be very difficult to reach the absolute values of neu¬ 
trino masses in terrestrial experiments. Future cosmological measurements could be a 
possibility in this case (see Ref. 2081 1. although the inferred information is somewhat 
model dependent. If m± is much smaller than ■>/Am 2 tm ~ 5 • 10~ 2 eV and the neutrino 
mass spectrum is characterized by an inverted hierarchy, information about the abso¬ 
lute values of neutrino masses can be obtained from the future experiments searching for 
neutrinoless double-/? decay. If mi is much larger than a/A m 2 tm (practically degener- 
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ate neutrino mass spectrum), there is a possibility that m\ will be measured in future 
/3 -decay experiments. Of course, progress in experimental techniques could bring some 
unexpected surprises for the solution of the exciting problem of the absolute neutrino 
masses. 


10 Note added in proof 

After completion and submission of this review, the Wilkinson Microwave Anisotropy 
Probe (WMAP) Science Team issued a most impressive set of data on the cosmic tem¬ 
perature angular power spectrum and the temperature-polarization angular power spec¬ 
trum [332]. Analysis of the data by the WMAP team has strengthened the notion that 
a Standard Model in Cosmology is emerging. Indeed, a flat universe filled with vacuum 
energy, dark matter, baryons and structure arisen from a nearly scale-invariant spectrum 
of primordial fluctuations seems to fit not only the WMAP high precision data but also 
smaller scale CMB data, large scale structure data, and Supernova la data. Furthermore, 
it is fully consistent with a much wider set of astronomical data such as the baryon to 
photon ratio derived from observations of D/H in distant quasars, the Hubble parameter 
measurement by the HST Key Project, the size of mass fluctuations obtained from galaxy 
clusters studies, and the inferred ages of stars |340j . 

As was stated in Section 6, neutrino mass has an impact on cosmological and astro- 
physical observables, notably on the suppression of power at low scales of the matter 
density perturbation spectrum (see Section 6.3). This is in fact the key ingredient in 
the determination of neutrino mass bounds from cosmology, as neutrinos with masses of 
the order of 1 eV are easily mistaken at recombination for cold dark matter and thus 
have small influence on the CMB. Hence the new WMAP data [339.1340) permits an 
improvement of these bounds only indirectly, as the parameters that define the overall 
cosmological picture become more precise. So, in order to constrain neutrino masses be¬ 
yond what was reached in previous analysis (see Section 6.5) the WMAP team undertook 
a study of their data combined with other CMB data that explore smaller angular scales 
(ACBAR |341| and CBI [ 168j ) and, most importantly, with the large scale structure 2dF- 
GRS data, and with the power spectrum recovered from Lyman a forest measurements. 
As a result of their analysis they find: 

0.71 eV (95%C.L.), (10.1) 

i 

which implies, for 3 degenerate neutrino species, rrii < 0.23 eV . This is about a factor 
three better than the bounds ( 6 . 8 ) and (6.9) (which, moreover, had to assume strong 
priors on the Hubble constant and the matter density). The precise role played by 
priors in lifting degeneracies among cosmological parameters when using WMAP data to 
constrain neutrino masses has been analyzed in y>42i and |343] . Two priors turn out to be 
particularly important: the prior on the matter density Q m (not surprisingly, because the 
suppression of growth of matter fluctuations at small scales depends on the ratio Q u /Q m ) 
and the prior on the Hubble parameter h. Hence, the WMAP results (where only flat 
models and h > 0.5 were allowed) Q m h 2 = 0.14 ± 0.02 and h = 0.72 ± 0.05 |340j are 
crucial in setting the neutrino mass limit given above. In their study, the authors in [342| 
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include an analysis of the 2dFGRS data complemented with simple gaussian priors from 
WMAP piD i: that leads to the following bound: 

J^rrii<l.leV (95%C.L.). (10.2) 

i 

In [ 343 ] Hannestad performs a search for neutrino mass bounds in three stages. In 
his most conservative analysis, where only WMAP data and 2dFGRS are being used, he 
obtains 

J2 m i< 2.12eV (95% C.L.). (10.3) 

i 

In a second stage, he includes extra CMB data from the compilation in [ 344] that incor¬ 
porates data at large Vs. He finds the considerably more restrictive limit: 

< 1.20eV (95% C.L.). (10.4) 

i 

The reason for improvement can be traced to the fact that the scales probed by the 
CMB data at high Vs [3441 overlap significantly with the scales under scrutiny in the 
2dFGRS survey. As a consequence a normalization of the 2dFGRS power spectrum from 
the CMB data is possible. In the final analysis Hannestad adds the constraints on h from 
the Hubble Space Telescope key project HST D23 and the constraints on matter density 
from the SNIa project GEt As a result, the bound on the total neutrino mass is, 

Y m % < 1.01 eV (95% C.L.). (10.5) 

i 

We see from these numbers that neither _ 342j nor [ 5131 can match the limit obtained 
by the WMAP team, a fact that the authors in these latter studies attribute mainly to not 
using Lycc data whereas the WMAP team does. Given that the extraction of the matter 
power spectrum from the Lyo forest involves complex numerical simulations, both [342 ] 
and .‘33| argue that one is on a safer position if these data are not used. 
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